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=  see  (11-27 ) ,  ( 11-28),  (11-29)  and  (11-30) 

=  gain  value  of  attitude  control  channel 
=  gain  value  of  control  damping 
=  see  (11-27),  (11-28),  (11-29)  and  (11-30) 

=  gain  value  of  flow  direction  channel 
=  see  (11-27),  (11-28),  (11-29)  and  (11-30) 
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=  base  diameter  of  missile 
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=  acceleration  of  gravity 
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=  instantaneous  immersion 

=  mass  moment  of  inertia  of  missile  about  its  c.g. 
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p  =  air  density 

t  =  time 

=  amplitude  of  first  sloshing  mass  relative  to  tank  wall 

§  =  amplitude  of  second  sloshing  mass  relative  to  tank  wall 

0  -  pitch  angle 

0.  =  indicated  attitude  angle 
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Subscripts 
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( 1 )  =  denotes  differentiation  with  respect  to  x 
( * )  =  denotes  differentiation  with  respect  to  time 


ABSTRACT 


,  The  wind  induced  responses  of  the  Saturn  C-5  without  fins  are  cal¬ 
culated  with  three  aerodynamic  representations .  The  most  accurate  representa¬ 
tion  uses  unsteady  aerodynamics  and  accounts  for  penetration  into  the  gusts . 

The  second  uses  pseudo  steady  aerodynamics  and  accounts  for  penetration.  The 
third  uses  pseudo  steady  aerodynamics  and  assumes  equal  wind  cross  flows  over 
the  missile  length. 

The  responses,  which  include  two  sloshing  and  two  bending  modes,  are 
affected  by  penetration  to  a  detectable  but  insignificant  degree .  The  use  of 
unsteady  aerodynamics  causes  very  little  change .  The  conventional,  third  aero¬ 
dynamic  representation  provided  in  these  cases  slightly  conservative  (large) 
responses . 
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SUMMARY 


f Missile  responses  to  winds  are  usually  calculated  using  pseudo  steady, 
slender  body  aerodynamics  and  the  assumption  that  wind  induced  crossflows  are 
equal  at  all  stations  along  the  missile.  However,  it  has  been  recognized  that 
the  time  delays  associated  with  missile  penetration  into  a  gust  could  signifi¬ 
cantly  affect  the  actual  forcing  functiohS,  especially  those  for  the  bending 
modes . 

The  wind  induced  responses  of  the  Saturn  C-5  (without  fins)  are  cal¬ 
culated  and  compared  for  three  aerodynamic  representations.  The  most  accurate 
representation  accounts  for  the  time  delays  of  penetration  and  uses  unsteady 
aerodynamics  so  that  growth  of  lift  smoothing  and  delays  are  included.  The 
second  representation,  called  pure  penetration,  uses  pseudo  steady  aerodynamics 
but  includes  penetration  delays.  The  third  representation  is  the  conventional 
one,  defined  here  as  instantaneous  immersion. 


The  analytical  model  includes  t^e v/f olljowing  degrees  of  freedom: 
translation,  rotation,  first  and  second  betiding"/' xwo  sloshing  modes,  and  con¬ 
trol.  Frozen  coefficients  are  used  with  each  set  applicable  to  a  10-sec. 
flight  time. 


Responses  are  calculated  for  the  unit  step  and  unit  impulse  gusts. 
These  responses  are  used  in  a  Duhamel  integration  to  obtain  the  responses  to 
two  wind  profiles  defined  by  values  at  25-meter  altitude  increments. 

The  incorporation  of  penetration  effects  causes  a  detectable  but 
insignificant  change  in  the  responses  calculated  for  the  Saturn  C-5.  The  ad¬ 
dition  of  growth  of  lift  causes  a  very  small  change.  The  conventional,  in¬ 
stantaneous  immersion,  responses  are  conservative  (large)  for  the  cases  com¬ 
puted. 


The  present  calculations  do  not  provide  an  example  of  significant 
response  changes  due  to  penetration  and  lift  growth.  However,  it  appears  that 
the  possibility  and  nature  of  significant  effects  could  be  detected  in  a  com¬ 
parison  of  the  unit  step  and  impulse  responses  based  on  instantaneous  immer¬ 
sion  and  pure  penetration. 

It  is  recommended  that  the  response  comparisons  be  extended  to  the 
third  and  fourth  bending  modes  and  to  a  Saturn  C-5  with  (simulated)  fins.  I 
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I*  INTRODUCTION 


The  calculation  of  missile  responses  due  to  winds  requires  the  in¬ 
corporation  of  a  satisfactory  description  of  the  aerodynamic  forces  into  the 
response  calculations.  In  a  previous  program  [l] ,  two  aerodynamic  effects  were 
pointed  out  which  have  been  neglected  in  previous  investigations  of  missile 
response  to  horizontal  winds  or  gusts.  The  first  is  the  gust  penetration  effect 
and  the  second  is  the  lift  growth  lag  due  to  aerodynamic  inertia.  In  the  same 
reference  these  effects  were  incorporated  in  rigid  body  indicial  and  impulsive 
aerodynamic  loads  based  on  slender  body  theory.  Some  numerical  results  were 
presented  for  a  multistage  missile  entering  a  unit  step  and  impulse  gust.  The 
results,  however,  were  not  utilized  in  response  calculations. 

Some  effects  of  penetration  and  growth  of  lift  were  postulated  in  [l ]  . 
For  rigid  body  motions,  the  effect  will  be  to  filter  out  the  high  wind  fre¬ 
quencies.  On  the  other  hand,  when  missile  bending  modes  are  included,  there 
is  the  possibility  of  augmenting  some  wind  frequency  responses. 

The  project  reported  here  is  an  investigation  of  the  effects  and 
importance  of  penetration  and  growth  of  lift  on  missile  responses  when  both 
rigid  body  and  bending  modes  are  considered.  The  Saturn  C-5  configuration, 
excluding  fins,  is  used  in  the  study. 


II.  BASIC  GOAL  OF  THE  INVESTIGATION 


The  basic  goal  of  the  work  presented  in  this  report  was  to  determine 
the  effects  and  importance  of  penetration  and  growth  of  lift  on  missile  re¬ 
ponses.  To  accomplish  this  goal,  the  investigation  was  divided  into  two 
phases: 


1.  The  indicial*  and  impulsive**  responses  of  the  vehicle  were  com¬ 
puted  using  three  aerodynamic  considerations:  (l)  penetration  with  lift  growth 
effects,  (2)  pure  penetration,  and  (3)  instantaneous  immersion.  Slender  body 
theory  was  used  throughout  the  investigation. 

2.  The  wind- induced  responses  of  the  vehicle  for  the  three  aero¬ 
dynamic  environments  were  computed  from  the  impulsive  responses  using  the 
Duhamel  superposition  integral. 


*  Response  to  a  unit  step. 

**  Response  to  a  unit  impulse . 
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In  both  phases,  numerical  comparisons  were  made  of  responses  calcu¬ 
lated  both  with  and  without  the  effect  of  penetration  and  lift  growth.  A 
mathematical  model  of  the  Saturn  C-5  configuration,  excluding  fins,  was  used 
in  the  analysis. 


III.  DISCUSSION  OF  THE  INVESTIGATION 


A.  Description  of  Saturn  C-5  Model  Used 


The  Saturn  C-5  configuration  was  used  in  the  investigation  of  the 
effects  and  importance  of  penetration  and  growth  of  lift  on  missile  response. 
The  equations  of  motion  of  the  missile  system  are  presented  in  Appendix  I. 
Seven  generalized  coordinates  were  considered  in  the  response  calculations; 
lateral  translation,  yQ  ;  rotation,  0  ;  the  first  two  bending  coordinates, 

%  and  7)2  ;  two  sloshing  coordinates,  ^  an^  ?2* **  >  and  con'trol  deflection, 
.  The  control  system  considered  in  this  report  utilizes  both  an  attitude 
reference  control  and  a  flow  direction  indicator. 

The  equations  are  valid  for  a  swivel  engine  controlled  vehicle  where 
the  swiveled  engines  account  for  four-fifths  of  the  total  thrust  force.  The 
missile  and  atmospheric  parameters  appearing  in  the  equations  are  considered 
constant  in  predetermined  time  or  altitude  intervals. 

The  actual  Saturn  C-5  configuration  contains  engine  shrouds  and  fins 
located  on  the  aft  section  of  the  vehicle.  The  effects  of  these  empennages 
were  neglected  in  the  analysis . 


B.  Description  of  Indicial  and  Impulsive  Aerodynamic  Forces  Used 

The  development  of  the  transient  and  quasi-steady,  generalized  aero¬ 
dynamic  forces  resulting  from  a  unit  step  and  impulsive  wind  profile  is  pre¬ 
sented  in  Appendix  II  of  this  report .  The  development  is  based  on  slender 
body  theory.  The  forces  corresponding  to  rigid  body  and  bending  coordinates 
are  presented  first  for  a  general  wind  profile.  The  indicial  and  impulsive 
forces  are  then  derived. 


*  The  two  sloshing  coordinates  considered  are  associated  with  the  fundamental 
fluid  motion  in  the  lox  and  fuel  tanks  located  in  the  booster  stage  of 
the  C-5. 

**  The  actual  engine  deflection,  0^  ,  and  the  control  deflection,  pc  ,  are 
considered  to  be  equal. 
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The  wind- induced  forcing  functions  which  are  compared  in  this  report 
have  both  a  geometric  and  an  aerodynamic  aspect. 

In  the  geometric  consideration  two  cases  are  used.  In  the  simplest 
case,  instantaneous  immersion,  all  stations  along  the  missile  are  assumed  to  be 
immersed  in  the  same  wind-induced  crossflow,  namely  the  wind  crossflow  occur¬ 
ring  at  the  nose.  The  more  accurate  geometric  representation,  called  penetra¬ 
tion,  assigns  to  each  station  along  the  missile  the  wind  crossflow  which  exists 
at  the  altitude  occupied  by  the  station.  With  penetration,  the  missile  nose 
enters  a  side  gust  first  and  in  subsequent  time  successive  stations  along  the 
missile  length  move  into  the  crossflow. 

Two  representations  of  the  aerodynamics  are  used;  they  are  quasi¬ 
steady  and  transient.  In  the  quasi-steady  representation  the  airforces  at  a 
missile  station  are  those  which  would  exist  if  the  local  crossflow  persisted 
unchanged  for  an  extended  time.  The  transient  representation  is  based  on  the 
theory  of  unsteady  motion  of  slender  bodies  and  includes  the  growth  of  lift 
with  time. 

Three  types  of  wind- induced  forcing  functions  are  assembled  using 
combinations  of  the  geometric  and  aerodynamic  representations .  The  simplest 
type  is  called  instantaneous  immersion  and  uses  the  instantaneous  immersion 
geometric  representation  with  quasi-steady  aerodynamics*  A  more  accurate  type, 
called  pure  penetration,  uses  penetration  geometries  and  quasi-steady  aero¬ 
dynamics  .  The  most  accurate  aerodynamic  forcing  functions  are  called  penetra¬ 
tion  with  lift  growth.  These  latter  functions  use  the  penetration  geometries 
with  transient  aerodynamics.* 

The  simpler  function  types,  instantaneous  immersion  and  pure  penetra¬ 
tion,  can  be  obtained  from  penetration  with  lift  growth  {see  Appendix  II  }  . 

The  development  of  the  indicial  and  impulsive  transient  and  quasi¬ 
steady  generalized  aerodynamic  forces  for  rigid  body  motion  follows  the  work 
presented  in  CG  and  Q0  . 

The  derivation  of  the  indicial  and  impulsive  transient  and  quasi¬ 
steady  aerodynamic  forces,  Q~  ,  corresponding  to  the  bending  coordinates. 


*  The  crossflows  induced  by  missile  motions  are  in  all  cases  treated  with 
quasi-steady  aerodynamics.  These  crossflows  are  small  in  comparison  to 
the  wind- induced  crossflows  and  appear  in  the  left  hand  side  of  the 
equations  of  motion  £see  Appendix  I  y  . 
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TlX(f  ,  requires  a  description  of  the  mode  shapes  of  the  missile.  A  considerable 
savings  in  computation  can  be  obtained  if  the  mode  shapes,  Ym( x )  ,  are  approxi¬ 
mated  by  polynomials.  The  components  of  corresponding  to  the  constant 

and  linear  terms  of  the  polynomials  can  then  be  rewritten  in  terms  of  the  wind- 
induced  aerodynamic  force  expressions  corresponding  to  the  rigid  body  coordi¬ 
nates.  In  addition,  the  polynomial  representations  for  Ym(x)  need  only  apply 
to  specific  regions  of  the  missile  length,  since  the  kernel  of  the  integrals 

describing  Qr,  take  on  values  only  over  the  conic  sections  of  the  missile. 

'm 

The  use  of  polynomial  approximations  for  Ym(x)  does  not  detract  from  a  general 
approach  to  numerical  solution,  since  the  important  mode  shapes  of  the  Saturn 
C-5  configuration  were  efficiently  described  by  low  order  polynomials  in  the 
regions  of  interest. 

Mode  shape  segments,  corresponding  to  the  conical  regions  of  the  C-5, 
were  fitted  with  quadratic  polynomials  of  the  form 

Ym(x)  =  Aim  +  BjLmx  +  ^  imx^ 

where  the  coefficients  A  ,  B  and  C* **  are  considered  constant  for  a  discrete 
flight  time  or  altitude. 

The  subscripts  m  and  i  designate  a  specific  mode  and  conical 
region,  respectively.  Four  conical  regions  (i  =  1,2, 3, 4)  were  considered  in 
the  analysis  {see  Appendix  II  }  .  For  the  C-5  configuration  ,  the  quadratic 
polynomial  mode  shapes  yielded  values  which  are  within  1  per  cent  of  the  actual 
mode  deflections . 


*  Penetration  and  lift  growth  lag  effects  are  potentially  important  when 

missile  bending  modes  are  considered  in  an  analysis  of  missile  response. 
Since  the  wind-induced  aerodynamic  forces  corresponding  to  the  bending 
coordinates  were  not  developed  in  {lj ,  it  was  necessary  to  formulate  these 
forces  in  this  report.  These  forces  were  derived  using  slender  body  theory 
and  include  the  geometric  and  aerodynamic  aspects  discussed  above . 

**  numerical  values  of  these  coefficients  were  computed  for  the  first  four 
bending  modes  (m  =  1,2, 3,4)  between  flight  times  of  10  and  140  sec.  at 
10-sec.  intervals.  However,  values  for  only  the  first  two  bending  modes 
between  30  and  100  sec.,  inclusive,  were  used  in  the  numerical  computation 
of  the  indicial  and  impulsive  responses. 
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The  fundamental  expressions  required  to  compute  the  wind- induced 
forcing  functions  for  the  Saturn  C-5  are  given  in  Appendix  II.  The  equation 
numbers  of  the  indicial  and  impulsive  forcing  functions  are  given  below  for 
each  of  the  three  aerodynamic  considerations. 


Aerodynamic  Consideration  Equation  Numbers 

1.  Penetration  with  lift  growth 

a.  Indicial 

b .  Impulsive 

2 .  Pure  penetration 

a.  Indicial 

b .  Impulsive 

3.  Instantaneous  immersion 

a.  Indicial  ( 11-66),  ( 11-67),  (11-68) 

b.  Impulsive  ( 11-69),  (11-70),  (11-71) 


(11-13),  (11-20),  (11-31) 
(11-38),  (11-39),  (11-40) 

(11-51),  (11-53),  (11-55) 
(11-57),  (11-59),  (11-61) 


For  purposes  of  general  interest,  plots  of  the  C-5  normal  force, 
first  bending  moment  and  second  bending  moment  for  a  unit  step  and  impulse 
side  wind  are  given  in  this  section.  The  curves  are  presented  for  a  Mach 
number  of  1.345  (70~sec.  flight  time).  The  consideration  of  penetration  ef¬ 
fects  causes  the  generalized  forces  to  be  distributed  over  time.  These  time- 
dependent  forcing  functions  are  called  growth  functions.  The  indicial  normal 
force ,  first  bending  moment  and  second  bending  moment  growth  functions  are 
given  in  Figs.  1,  2  and  3,  respectively.  The  impulsive  normal  force,  first 
bending  moment  and  second  bending  moment  growth  functions  are  given  in  Figs . 

4,  5  and  6,  respectively.  The  effects  of  penetration  with  lift  growth,  pure 
penetration  and  instantaneous  immersion  on  the  generalized  forces  can  be 
easily  seen  from  the  figures.  The  impulsive  forces  resulting  from  the  effects 
of  instantaneous  immersion  are  not  presented  since  they  can  only  be  described 
from  a  limit  consideration. 
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-  Indicial  2nd  Bending  Moment  Growth  Functions  vs.  Time  in  Seconds 


h 

M 

[  _ Pure  Penetration 

'  — _ Penetration  with  Lift  Growth 


4  -  Impulsive  Normal  Force  Growth  Functions  vs.  Time  in  Seconds 


rulsive  Second  Bending  Moment  Grovth  Function  xlO~^  ^  Impulsive  First  Bending  Moment  Grovth  Function  xlO 


_  Pure  Penetration 

-  Penetration  with  Lift  Grovth 

M  -  1.345 


0.02  0.04  0.06  0.08  0.12  0.14  0.16  0.18  0.20  0.22  0.24  0.26  0.28  0.30 


Impulsive.  1st  Bending  Moment  Grovth  Functions  vs.  Time  in  Seconds 


0.02  0.04 


—  Pure  Penetration 

-  Penetration  vith  Lift  Grovth 

M  =  1.345 


0.06  0.08  0.10  0.12  0.14  0.16  0.18  0.20  0.22  0.24  0.26  0.28  0.30 

Time ,  t 


Impulsive  2nd  Bending  Moment  Grovth  Functions  vs.  Time  in  Seconds 


The  half  cross  section  of  the  C-5  configuration  is  given  at  the  bottom 
of  each  of  the  figures  for  convenience  of  interpretation  of  the  growth  functions. 
At  t  =  0  ,  the  first  region  (escape  tower)  penetrates  the  gust  front  and  the 
generalized  force  buildup  begins.  The  times  t  —  0.0216,  0.0500,  and  0.0930 
sec.  correspond  to  the  successive  gust  encounters  of  the  second,  third  and 
fourth  conic  regions.  The  pure  penetration  and  penetration  with  lift  growth 
curves  reached  steady-state  conditions  at  r  =  0.1065  sec.  and  t  =  0.2695  sec., 
respectively.  The  time  for  total  immersion  of  the  vehicle  was  0.2629  sec. 

For  reasons  of  brevity,  a  detailed  description  of  the  growth  functions 
will  not  be  given  in  this  report.  The  reader  is  referred  to  [1]  for  some  gen¬ 
eral  comments  concerning  the  rigid  body  growth  functions.  The  general  form  of 
the  bending  moment  growth  functions  are  very  similar  to  those  of  the  normal 
force  functions.  The  bending  moment  growth  functions,  however,  reflect  the 
character  of  the  mode  shapes.  The  negative  growth  functions  for  the  fourth 
region  in  Fig.  5  and  the  first  and  second  regions  in  Fig.  6  correspond  to  re¬ 
gions  of  negative  displacements  of  the  first  and  second  bending  modes,  respec¬ 
tively.  Figures  2  and  3  reflect  an  integrated  effect  of  the  growth  functions 
in  Figs.  5  and  6,  respectively. 


C.  Method  of  Solution  for  Indicial  and  Impulsive  Responses  of  Saturn  C-5 

The  method  of  Runge-Kutta  was  used  to  calculate  the  impulsive  and 
indicial  responses  of  the  Saturn  C-5.  A  fourth  order*  Runge-Kutta  integration 
equation  was  used  to  numerically  integrate  the  equations  of  motion.  Saturn 
C-5  missile  response  calculations  were  obtained  for  six  different  sets  of  aero¬ 
dynamic  environments.  These  six  sets  of  forcing  functions  {see  Appendix  III, 
Table  VII  }  contain  the  indicial  and  impulsive  aerodynamic  forces  resulting 
from 


1.  Penetration  with  lift  growth  effects, 

2.  Pure  penetration  effects,  and 

3.  Instantaneous-immersion  effects. 


*  The  Runge-Kutta  integration  expression  used  is  considered  to  be  of  fourth- 
order  accuracy  {  see  C3]  }  . 
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The  Runge-Kutta  integration  process  requires  knowledge  of  the  initial 
conditions  (or  starting  values).  Since  the  system  of  equations  used  in  the 
investigation  can  be  written  as  a  second  order  set,  ^see  (il-ll)  through 
( 11-19)  only  the  initial  conditions  for  the  dependent  variables  and  their 

first  derivatives  are  needed.  The  initial  conditions  for  each  of  the  six  sets 
of  forcing  functions  were  found  through  use  of  Laplace  transform  techniques 
«£see  Appendix  III j  . 

The  indicial  and  impulsive  responses  of  the  C-5  were  computed  at  dis¬ 
crete  flight  times  or  altitudes  from  a  system  of  linear,  differential  equations 
with  constant  coefficients*  {see  (III-1)}.  Each  set  of  constant  coefficients, 
and  consequently  the  associated  responses,  is  applicable  in  a  specified  time 
or  altitude  interval.**  A  computer  program  used  to  obtain  the  response  calcu¬ 
lations  is  discussed  in  Vol.  III. 

In  each  flight  time  or  altitude  band,  the  responses  were  computed 
for  a  real  time  interval  of  20  sec .  Extremely  fine  increments  of  response  time 
were  used  in  computing  the  forcing  functions  and  the  indicial  and  impulsive 
responses  of  the  deflections  and  their  first  derivatives .  Coarser  increments 
were  used  in  the  printing  and  plotting  of  the  output.  Only  the  deflections 
were  plotted  on  the  SC -4020. 

Twenty- second  response  records  were  computed  to  insure  that  the  re¬ 
sponses  had  achieved  satisfactory  steady-state  values.  For  the  flight  times 
considered,  all  of  the  impulsive***  responses  had  converged  to  1  per  cent  or 
less  of  their  maximum  values  at  the  end  of  20  sec . 

The  calculation  of  20-sec.  response  records  does  not  violate  the  as¬ 
sumption  of  constant  coefficients  as  might  be  expected.  From  a  subsidiary  re¬ 
sponse  analysis,  it  was  found  that  if  the  sloshing  degrees  of  freedom  were 
removed  from  the  system,  the  responses  of  the  remaining  system  (excluding 
translation)  would  achieve  satisfactory  steady-state  values  within  7  or  8  sec. 


*  Previous  work  has  shown  that  it  is  permissible  to  use  sets  of  constant 
coefficients . 

**  Response  calculations  were  obtained  at  every  10  sec.  of  flight  time  between 
30  and  100  sec.  The  applicable  interval  of  flight  time  for  each  set  of 
these  responses  is  taken  as  10  sec.  Thus,  for  example,  the  30-sec.  re¬ 
sponse  calculations  (obtained  using  the  coefficients  at  30  sec.)  are  con¬ 
sidered  valid  between  flight  times  of  25  and  35  sec. 

***The  decay  of  the  indicial  responses  will  be  discussed  later  in  this  sec¬ 
tion. 
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of  response  time.  The  sloshing  modes  of  the  C-5  are  very  slightly  damped  in 
comparison  to  the  other  modes  of  the  system  (except  translation).  Therefore, 
since  we  are  considering  a  coupled  system,  the  response  calculations  for  rota¬ 
tion,  first  and  second  bending  and  control  after  7  or  8  sec.  reflect  the  in¬ 
fluence  of  the  sloshing  modes.  This  influence  is  basically  governed  by  that 
part  of  the  equation  of  motion  which  contains  time  invariant  coefficients. 

The  selection  of  the  increment  size  used  in  the  numerical  calculation 
of  the  responses  was  dictated  by  two  requirements: 

1.  The  increment  size  should  be  sufficiently  small  as  to  permit  an 
accurate  calculation  of  the  aerodynamic  forces. 

2.  The  increment  size  should  permit  32  calculated  response  values 
per  cycle  for  the  highest  frequency  component  of  the  system. 

Based  on  the  above  requirements,  the  increment  sizes.  At  ,  used  in 
the  response  calculations  are  presented  in  Table  I. 

TABLE  I 

TABULATION  OF  INCREMENT  SIZES  USED  IN  INDICIAL 
AND  IMPULSIVE  RESPONSE  CALCULATIONS- 


Flight  Time 
(sec . ) 

AT1 

(sec .  ) 

Ti 

(sec . ) 

At2 

( sec . ) 

T2 

(sec .  ) 

At3 
(sec .  ) 

30 

0.0005 

0.02 

0.001 

0.75 

0.004 

40 

0.0005 

0.02 

0.001 

0.60 

0.004 

50 

0  .0005 

0.02 

0.001 

0.50 

0.004 

60 

0.0005 

0.02 

0.001 

0.45 

0.003 

70 

0.0005 

0.03 

0.001 

0.40 

0.003 

80 

0 . 0005 

0.03 

0.001 

0.40 

0.004 

90 

0.0005 

0.04 

0.001 

0.35 

0.004 

100 

0.0005 

0.05 

0.001 

0.35 

0.004 
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The  response  times ,  t  ,  at  which  the  increment  size  was  changed*  are 
also  given,  At^  was  used  in  the  calculations  from  0  to  ;  At2  was  used 
from  ti  to  T2  ;  and  AT3  was  used  from  t2  to  the  end  of  the  response 
calculations . 

Extremely  small  increments  (  Arp  =  0,0005  sec.)  are  used  at  the  begin¬ 
ning  of  the  response  calculations  to  adequately  describe  the  high-frequency 
components  of  the  control  system.  The  Tj_  values  correspond  to  the  times  at 
which  these  frequency  components  are  negligible.  The  T2  values  correspond  to 
conservative  estimates  of  the  times  when  the  aerodynamic  forcing  functions 
reach  steady-state  conditions. 

Since  a  voluminous  amount  of  indicial  and  impulsive  response  data 
were  generated,  only  a  representative  quantity  of  these  data  is  presented.  Plots 
of  the  indicial  and  impulsive  responses  of  first  bending,  7^  ,  second  bending, 
T]2  >  and  control  deflection,  $c  (noted  by  first  control)  are  given  in  Figs.  7 
through  24  for  a  flight  time  of  70  sec.  The  impulsive  responses  are  presented 
in  Figs.  7  through  15;  the  indicial  responses  are  given  in  Figs.  16  through  24. 
For  each  coordinate,  the  first  figure  reflects  the  effect  of  instantaneous  im¬ 
mersion;  the  second  figure  reflects  the  effect  of  pure  penetration;  and  the 
third  reflects  the  effect  of  penetration  with  lift  growth. 

The  numerics  presented  alongside  the  plots  in  each  figure  pertain  to 
an  analysis**  of  the  zero  crossings,  maximum  and  minimum  values  of  the  respec¬ 
tive  coordinate  response.  The  format  of  the  information  presented  is  as  fol¬ 
lows  : 


Response 

Time 


Response  Response 

Value  Maximum,  Minimum  or 

Crossing 


The  maximum,  minimum  and  zero  crossing  values  are  the  calculated  points  which 
precede  (in  time)  the  event. 


*  It  was  found  expedient,  from  the  standpoint  of  conserving  computer  running 
time,  to  increase  the  increment  size,  when  permissible,  during  the 
response  calculations. 

**  The  write-up  on  the  computer  program  used  to  analyze  the  indicial  and 
impulsive  responses  is  presented  in  Vol.  III. 
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Impulsive  Response  of  1st  Bending  Considering  Instantaneous  Immersion 
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Fig.  8  -  Impulsive  Response  of  1st  Bending  Considering  Pure  Penetration 
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Fig.  10  -  Impulsive  Response  of  2nd  Bending  Considering  Instantaneous  Immersion 
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Fig.  12  -  Impulsive  Response  of  2nd  Bending  Considering  Penetration  with  Lift  Growth 
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Fig.  14  -  Impulsive  Response  of  1st  Control  Considering  Pure  Penetration 
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Fig.  15  -  Impulsive  Response  of  1st  Control  Considering  Penetration  with  Lift  Growth 
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Fig.  16  -  Indicial  Response  of  1st  Bending  Considering  Instantaneous  Immersion 


25 


Fig.  18  -  Indicial  Response  of  1st  Bending  Considering  Penetration  with  Lift  Growth 
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Fig.  19  -  Indie ial  Response  of  2nd  Bending  Considering  Instantaneous  Immersion 
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Fig.  21  -  Indicial  Response  of  2nd  Bending  Considering  Penetration  with  Lift  Growth 
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Fig.  22  -  Indicial  Response  of  1st  Control  Considering  Instantaneous  Immersion 
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Fig.  23  -  Indicial  Response  of  1st  Control  Considering  Pure  Penetration 
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Fig.  24  -  Indicial  Response  of  1st  Control  Considering  Penetration  with  Lift  Growth 


The  discussion  of  the  effects  and  importance  of  penetration  and  growth 
of  lift  on  the  indicial  and  impulsive  responses  of  the  vehicle  will  be  presented 
in  a  later  section. 

There  is,  however,  some  information  readily  derivable  from  the  indicial 
and  impulsive  responses  which  has  a  direct  bearing  on  the  calculation  of  the 
wind-induced  responses.  The  evaluation  of  the  wind-induced  responses  can  be 
determined  from  either  the  indicial  or  impulsive  responses .  An  important  fea¬ 
ture  in  this  choice  is  the  convergence  characteristic  of  the  response.  From 
Figs.  1  through  18,  it  is  seen  that  the  indicial  responses  do  not  converge 
satisfactorily  within  the  20-sec.  record.  The  path  or  drift  root  component  has 
a  predominant  influence  on  the  slow  convergence  of  the  indicial  responses.  Hie 
impulsive  responses,  however,  are  not  as  greatly  influenced  by  this  component 
and  do  converge  satisfactorily  within  the  record.  Thus,  the  impulsive  responses 
were  selected  for  use  in  the  calculation  of  the  wind-induced  responses .  The 
evaluation  of  the  wind-induced  responses  is  discussed  in  the  next  section. 


D.  Wind-Induced  Responses 


This  section  presents  the  method  used  and  procedures  followed  in  com¬ 
puting  the  wind-induced  responses  of  the  Saturn  C-5  configuration.  A  sample  of 
the  numerical  calculations  is  presented  for  illustrative  purposes . 

The  responses  of  the  vehicle  to  any  side  wind  input  are  formulated  as 
Duhamel  integrals  £see  [4]  ]►  .  The  kernel  function  of  these  integrals  can  be 
written  in  terms  of  either  the  indicial  or  impulsive  responses  of  the  system. 
From  the  discussion  given  in  the  previous  section,  it  is  advantageous  to  com¬ 
pute  the  wind-induced  responses  in  terms  of  the  impulsive  responses. 

Thus,  from  (jT),  the  wind-induced  responses  are  given  by 


R(t)  =  J  Ri(t-T)vy(T)dT 

where  H1(t)  is  the  Ith  coordinate  response  to  a  unit  impulse  wind,  vy(T)  is 
an  arbitrary  side  wind  profile  and  R(t)  is  the  i^h  coordinate  response  to  the 
wind,  Vy(t). 
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The  wind  velocity  used  in  the  above  integral  was  available  as  a 
function  of  altitude.  The  profile  was  interpreted  in  terms  of  time  by  using  a 
constant  missile  velocity  associated  with  the  midpoint  of  the  altitude  band. 
This,  in  essence,  produced  a  slight  shift  in  the  effective  wind  frequencies  at 
the  extremes  of  the  altitude  interval. 

The  responses  of  the  Saturn  C-5  to  two  wind  profiles*  were  computed. 
The  wind  profiles  used  were:  the  West  -  East  components  of  8116  and  4768.** 

(in  this  report  8116  is  referred  to  as  profile  1;  and  4768  is  referred  to  as 
profile  2.)  These  profiles  had  been  measured  by  a  modified  spherical  balloon 
(Jimsphere)  -  radar  technique  and  were  tabulated  at  25-meter  increments.  The 
8116  profile  represents  a  moderately  severe  wind  environment  (maximum  wind 
velocity  of'  46.3  meters/sec  at  an  altitude  of  13,200  meters).  The  4768  profile 
represents  a  mild  wind  environment  (maximum  wind  velocity  of  17.8  meters/sec 
at  an  altitude  of  19,000  meters). 

The  vehicle  responses  to  the  two  profiles  were  computed  at  numerous 
altitude  points  during  the  flight  trajectory.  The  responses  of  the  following 
generalized  coordinates  were  computed: 

Translation,  yQ 

Rotation,  0 

First  bending, 

Second  bending, 

First  sloshing, 

Second  sloshing, 

Control  deflection, 


*  The  wind  data  were  supplied  by  George  C.  Marshall  Space  Flight  Center. 

**  Wind  values  for  profile  8116  were  given  between  altitudes  of  2,350  and 

19,525  meters;  values  for  profile  4768  were  given  between  450  and  19,700 
meters . 
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The  wind-induced  responses  were  evaluated  in  five  flight  time  hands. 
Each  hand  was  initially  chosen  to  encompass  10  sec.  of  flight  time  with  missile 
parameters  evaluated  at  the  hand  midpoint.  Subsequently,  each  of  the  actual 
calculations  was  extended  to  a  length  of  15  sec .  so  that  5  sec .  of  overlap  were 
available  for  comparison.  The  comparisons  were  a  practical  test  of  the  use  of 
frozen  coefficients.  The  agreement  of  response  values  in  the  overlap  regions 
was,  in  general,  satisfactory.  As  expected,  the  high-frequency  responses  showed 
better  agreement  than  low-frequency  responses . 

The  correspondence  between  the  altitudes  and  flight  time  bands  used 
in  the  numerical  calculation  of  the  wind-induced  responses  is  presented  in  Table 
II.  The  increment  used  in  each  band  to  calculate  the  responses  is  also  given. 
These  increments  were  chosen  to  provide  sufficient  calculated  points  to  ade¬ 
quately  define  the  highest  frequency  component. 


TABLE  II 

ALTITUDE  BANDS  USED  IN  CALCULATING  WIND-INDUCED  RESPONSES 


Bands 


Flight  Time 
( sec . ) 

Flight  Time 
Interval 
(sec . ) 

Altitude  Interval—/ 
(meters ) 

Response  Increment 
(meters ) 

50 

45  to  60 

3,594  to  7,104 

8.33 

60 

55  to  70 

5,787  to  10,200 

12.5 

70 

65  to  80 

8,573  to  13,940 

12.5 

80 

75  to  90 

11,987  to  18,361 

12.5 

90 

85  to  92.5 

16,063  to  19,525 

12.5 

a/  It  was  not  possible  to  obtain  an  adequate  description  of  the  slowly  con¬ 
verging  responses  at  the  lowest  end  of  some  of  the  intervals  because  the 
wind  profile  did  not  extend  to  a  sufficiently  low  altitude . 
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To  avoid  burdening  the  report  with  an  excessive  amount  of  data,  only 
a  representative  quantity  of  the  wind-induced  response  data  is  presented  at 
this  time.  Plots  of  the  first  bending,  %  ,  second  bending,  1)2  ,  and  control 
deflection,  pc  (noted  by  first  control)  responses  to  profile  8116  (noted  by 
profile  l)  are  given  in  Figs.  25  through  33  for  the  altitude  band,  8,573  to 
13,940  meters.  The  70-sec.  impulsive  responses  were  used  in  computing  these 
responses.  For  each  coordinate,  the  first  figure  reflects  the  effect  of  in¬ 
stantaneous  immersion;  the  second  figure  reflects  the  effect  of  pure  penetra¬ 
tion;*  and  the  third  reflects  the  effect  of  penetration  with  lift  growth . 

The  numerics  presented  alongside  the  plots  in  each  figure  pertain 
to  an  analysis  of  the  response  curves.  The  format  of  the  information  presented 
is  as  follows: 


Response  Response  Response 

Altitude  Value  Maximum,  Minimum  or 

Crossing 

The  maximum  and  minimum  values,  indicated  by  a  check  mark,  define  the  approxi 
mate  envelope  of  the  extreme  excursions  of  the  response.  These  values  were 
used  in  the  analysis  of  the  wind-induced  responses  discussed  in  the  next 
section. 


IV.  METHOD  USED  IN  ANALYSIS  OF  EFFECTS  AM)  IMPORTANCE  OF  PENETRATION 
AM)  GROWTH  OF  LIFT  ON  MISSILE  RESPONSE 


The  analysis  of  the  indicial  and  impulsive  responses  is  largely 
qualitative . 

The  three  aerodynamic  representations  yield  similar  responses  to  the 
indicial  wind.  The  responses  to  the  unit  impulse  wind  are  also  very  much  alike. 
Samples  of  these  responses  are  shown  in  the  section,  "Method  of  Solution  for  the 
Indicial  and  Impulsive  Responses."  The  most  noticeable  differences  occur  during 
the  highly  oscillatory  response  shortly  after  immersion  in  the  gust. 


*  An  undesired  translation  of  altitude  was  made  in  the  calculation  of  the 
responses  to  profile  1  considering  pure  penetration  effects .  The 
altitudes  recorded  in  the  plots  and  tabulations  should  be  reduced  by 
100  meters . 
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Fig.  26  -  Wind-Induced  Response  of  1st  Bending  Considering  Pure 

Penetration  -  70  sec.  F.T.  Band 
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Fig.  27  -  Wind-Induced  Response  of  1st  Bending  Considering  Penetration 

With  Lift  Growth  -  70  sec.  F.T.  Band 
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Fig.  28  -  Wind-Induced  Response  of  2nd  Bending  Considering  Instantaneous 

Immersion  -  70  sec.  F.T.  Band 


PR0FILE  It  ALTITUDE  VS.  RESP0NSE-PURE  PENETRAT I0N  70  SEC 

ETA2,  2ND  BENDING 
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Fig.  29  -  Wind-Induced  Response  of  2nd  Bending  Considering  Pure 

Penetration  -  70  sec.  F.T.  Band 


PR0FI LE  If  ALTITUDE  VS»  RES P 0 N  $JE -PEN ETRATI0N  W/  LIFT  GR0WTH  70  SEC 
E T A2 *  2ND  BENDING 
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Fig.  30  -  Wind-Induced  Response  of  2nd  Bending  Considering  Penetration 

With  Lift  Growth  -  70  sec.  F.T.  Band 
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Fig.  31  -  Wind- Induced  Response  of  1st  Control  Considering  Instantaneous 

Immersion  -  70  sec.  F.T.  Band 
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Fig.  33  -  Wind-Induced.  Response  of  1st  Control  Considering  Penetration 

With  Lift  Growth  -  70  sec.  F.T.  Band 


The  pure  penetration  and  penetration  with  lift  frowth  representations  yield 
slightly  smaller  oscillatory  excursions  than  instantaneous  immersion.  These 
differences  are  small  (not  more  than  10  per  cent  of  the  excursion  amplitude  or 
area)  and  it  is  difficult  to  predict  the  effect  on  calculated  responses  due  to 
actual  wind  profiles.  It  is  anticipated,  however,  that  the  aerodynamic  repre- 
senations  which  include  penetration  will  provide  smaller  oscillatory  excur¬ 
sions  .* 


A  second  feature  was  expected  and  verified  in  the  responses  to  impulse 
and  step  winds.  Those  responses  which  exhibit  a  long,  slow  approach  to  their 
steady-state  values  are,  during  this  approach,  nearly  independent  of  the  aero¬ 
dynamic  representation  used.  This  similarity  in  responses  indicates  that  the 
local  average  responses  to  wind  profiles  should  be  nearly  independent  of  the 
aerodynamic  representation  used. 

The  indicial  and  impulse  responses  are  seen  to  provide  clues  for  the 
analysis  of  wind  profile  responses.  In  the  wind-induced  responses  calculated 
with  the  three  aerodynamic  representations,  the  local  averages  are  expected  to 
be  similar  while  the  deviations  from  the  local  average  will  be  largest  using 
instantaneous  immersion.  These  ideas  about  the  wind  profile  induced  responses 
are  tested  using  formal  analytical  procedures. 

Responses  are  calculated  for  just  two  wind  profiles .  The  main  problem 
is  hew  to  use  this  small  amount  of  data  for  a  comparison  of  the  responses  from 
the  three  aerodynamic  representations.  The  comparisons  are  made  using  ex¬ 
tremals.  The  postulated  differences  in  excursions  from  the  local  average  are 
tested  by  comparing  the  distributions  of  extreme  excursions .  The  local  averages 
are  compared  using  the  average  of  the  envelope  of  points  which  in  pairs  define 
the  extreme  excursions . 

The  following  procedure  is  used  for  each  response  in  a  calculated 
altitude  interval.**  The  response  is  examined  and  those  extreme  points  are 
selected  which  appear  to  lie  on  or  near  envelope  curves  which  would  enclose 


*  The  reduction  in  excursions  is  important  because  maximum  responses  very 
often  will  be  due  to  excursions  or  deviations  from  a  "local  average." 

**  The  altitude  intervals  vary  from  3,510  to  6,375  meters.  Normally,  constant 
missile  characteristics  should  probably  not  be  used  for  intervals  which 
are  this  large.  However,  in  this  case,  the  additional  adjunct  wind  pro¬ 
file  can  be  considered  typical  of  the  altitude  region  and  thus  provides 
additional  typical  response  data. 
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the  response.*  The  values  of  the  response  at  these  points  are  read  and  then 
adjunct  values  are  differenced.  The  magnitudes  of  these  differences  are  ex^ 
treme  excursions  and  are  analyzed  as  a  sample  of  extremals  from  a  single 
population.  The  analysis  of  the  excursions  include  ordering,  transforming 
and  curve  fitting.  The  transformation  used  is 


where  y  is  called  the  reduced  variate,  §  =  -J5L  ,  the  accumulative  probabil- 

n+1 

ity  in  a  sample  of  n  maximums  and  m  is  the  order  from  the  smallest  (m  =  l). 

The  transformation  is  designed  to  provide  x  as  a  linear  function 
of  ym**  ,  where  x^  is  the  observed  ordered  maximum  in  the  sample  of  n 
maximums.  The  curve  fit  is  made  in  the  x,y  coordinates  using  the  method  of 
least  squares . 

A  local  average  response  is  calculated  as  the  unbiased  average  of  the 
extreme  point  values  initially  selected  from  the  response.  Vhen  an  odd  number 
of  points  is  selected,  the  number  of  maximums  and  minimums  are  unequal.  This 
bias  is  removed  with  a  weight  factor  of  one-half  on  the  first  and  last  points . 

The  analytical  techniques  described  here  constitute  unconventional 
employments  of  conventional  methods.  Ordinarily,  a  sample  of  extremes  is  ob¬ 
tained  by  taking  from  each  of  a  number  of  equal- sized  samples  the  largest 
(or  smallest)  value.  The  resulting  distribution  is  used  to  predict  probability 
of  occurrence  of  an  extreme  value  in  an  even  larger  number  of  similar-sized 
samples  taken  from  the  same  population.  An  associated  interpretation  provides 
the  return  period  which  is  the  expected  number  of  equal- sized  samples  which 
will  be  required  to  locate  an  assigned  maximum. 


*  Since  the  responses  being  compared  are  very  much  alike,  the  near  envelope 
points  selected  are  in  all  cases  equivalent  points. 

**  See  "Statistics  of  Extremes"  by  E .  J.  Gumbel,  Columbia  University  Press,  or 
"Statistical  Theory  of  Extreme  Values  and  Some  Practical  Applications"  by 
E.  J.  Gumbel,  National  Bureau  of  Standards,  Applied  Mathematics  Series  33. 
The  transformation  is  applicable  to  maximums  from  populations  with  proba¬ 
bility  functions  of  the  "exponential  type."  Many  important  distributions 
are  of  this  type  including  the  exponential  itself,  the  normal,  the  chi- 
square,  the  logistic,  and  the  log  normal. 
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In  the  present  analysis ,  the  basic  sample  size  is  the  altitude  in¬ 
crement  in  which  responses  are  calculated  with  frozen  missile  properties. 
However ,  not  just  one  maximum  excursion  is  taken  from  this  sample  of  the  calcu¬ 
lated  response.  The  number  of  excursions  taken  are  all  those  which  appear  to 
extend  to  envelope  curves.  This  selection  takes  the  attitude  that  the  response 
excursions  are  a  superposition  of  responses  involving  the  spectral  content  of 
the  wind  and  admittance  properties  of  the  vehicle.  With  this  viewpoint ,  the 
successive  "maximum  excursions"  may  be  considered  extremals  in  a  superposition 
of  independents ,*  thereby  justifying  the  use  of  the  analysis.  It  is  clear  that 
the  conventional  connotation  of  return  period  has  been  altered  here .  It  might 
be  possible  to  recover  this  predictive  capability  by  reinterpreting  sample 
size.  However,  this  predictive  use  of  the  present  analysis  is  not  recommended. 
Instead,  the  analysis  is  recommended  for  the  comparison  of  excursions  obtained 
using  the  three  aerodynamic  representations.  The  analysis  permits  the  compari¬ 
son  to  be  made  in  an  integrated  rather  than  point  fashion  and  uses  the  majority 
of  pertinent  data.  The  trend  with  y  ,  the  reduced  variate,  may  be  interpreted 
simply  as  the  nonlinear  dimension  of  increasing  numbers  of  samples  or  increas¬ 
ing  numbers  of  missile  transverses . 

The  calculation  of  an  average  from  the  average  of  the  extremes  is  not 
conventional  but  is  an  accepted  and  often  powerful  technique. 

Figure  34  shows  an  example  of  the  excursion  analysis  and  resulting 
curve  fits  plotted  on  paper  designed  for  this  purpose.  In  the  adjacent  Fig. 

35,  the  same  data  are  shown  on  rectilinear  coordinates,  the  observed  variate 
(ordinate)  and  the  reduced  variate  (abscissa)  which  is  the  transformed  cumula¬ 
tive  probability.  The  similar  curve  fits  for  the  analyzed  responses  are  shown 
on  rectangular  coordinates  in  Appendix  IV. 


*  The  problem  of  independence  in  the  primary  data  arises  in  the  analysis  of 
river  flows ,  There ,  it  is  an  inferred  assumption  that  the  flow  during 
each  24-hr.  period  is  an  independent  measure,  although  correlation 
between  the  flows  of  successive  days  is  easily  demonstrated.  The  theory 
of  the  extremes  is  still  very  successful  in  this  case. 
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V.  RESULTS  AND  DISCUSSION 


In  this  section  the  results  are  referenced  and  discussed  not  only 
for  their  immediate  import  but  also  as  a  guide  to  future  prediction  of  penetra¬ 
tion  and  lift  growth  effects. 

Samples  of  the  wind-induced  responses  are  shown  in  Figs.  25  through 
33  for  the  three  different  aerodynamic  representations.  The  most  important  re¬ 
sult  is  the  small  difference  between  the  responses  and  between  the  maximum  re¬ 
sponses.  The  same  result  is  observed  in  the  comparisons  for  other  responses  and 
flight  times.  Generally,  when  the  difference  in  maximum  responses  exceeds  1  per 
cent  the  instantaneous  immersion  results  are  conservative  (large).  The  differ¬ 
ences  in  maximum  responses  with  the  aerodynamic  representations  do  not  appear 
significant  in  an  engineering  sense.  For  instance,  with  profile  No.  1,  a  mod¬ 
erately  severe  profile,  instantaneous  immersion  provides  first  bending  which  is 
conservative  by  -0.7  per  cent  to  1.3  per  cent.  Second  bending  with  instantane¬ 
ous  immersion  is  conservative  by  - 1.0-  per  cent  to  2.3  per. cent,  while  the  first 
sloshing  response  is  conservative  by  -0.8  per  cent  to  4.8  per  cent. 

The  results  and  interpretations  which  follow  indicate  that  the  per 
cent  difference  may  be  larger  in  weak  wind  profiles  but  the  magnitude  of  the 
differences  will  remain  unin^ortant . 

The  results  indicate  that  the  responses  are  insignificantly  affected 
by  the  inclusion  of  penetration  effects  and  lift  growth  effects.  However,  the 
conclusion  can  be  drawn  only  for  the  missile  configuration  used  and  the  modes 
included.  The  possibilities  for  more  general  results  are  explored  next. 


Samples  of  the  indicial  and  impulsive  responses  are  shown  in  Figs.  7 
through  24.  There  is  a  small  but  noticeable  difference  between  the  responses 
which  do  and  do  not  include  penetration.  The  oscillatory  excursions  with  pene¬ 
tration  (PP  and  PWLG)  are  smaller  than  those  from  instantaneous  immersion.  This 
is  especially  true  for  the  bending  and  sloshing  modes.  It  was  then  anticipated 
that  the  same  type  of  difference  in  oscillatory  excursions  might  appear  in  the 
responses  to  wind  profiles.  An  analysis  of  the  wind-induced  excursions  was  de¬ 
signed  and  applied  as  described  in  Section  IV.  The  results  are  shown  in  Figs. 

42  through  89  (pp.  109  through  135)  where  the  anticipated  differences  are  ob¬ 
served.  In  these  figures  the  ordinate  is  the  magnitude  of  excursions  and  the 
abscissa  is  a  nonlinear  dimension  of  increasing  transverses  or  sample  sizes. 

Figures  42  through  89  illustrate  the  close  correspondence  between  re¬ 
sponses  with  pure  penetration  (PP)  and  penetration  with  lift  growth  (PWLG) .  The 
difference  between  responses  calculated  with  and  without  penetration  are  seen 
to  diminish  with  altitude.  This  would  be  expected  since  with  higher  speeds  the 
penetration  cases  are  approaching  instantaneous  immersion. 
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The  indicial  and  impulsive  responses  (Figs.  7  through  24)  also  indi¬ 
cate  an  important  similarity  in  responses  with  all  three  aerodynamic  represen¬ 
tations.  This  similarity  occurs  in  the  long  slow  approach  to  steady-state 
values.  The  local  average  values  of  the  wind-induced  responses  depend  largely 
on  this  long  tail  and  the  history  of  wind  inputs.  Thus,  local  average  responses 
to  winds  were  expected  to  be  about  equal  for  all  aerodynamic  representations. 
This  expectation  is  borne  out  by  the  averages  presented  in  Figs.  90  through  99 
(pp.  136  through  141). 

The  differences  and  similarities  in  the  indicial  and  impulsive  re¬ 
sponses  seem  to  have  their  logical  counterparts  in  the  wind-induced  responses 
taken  as  a  whole.  However,  since  main  interest  attaches  to  the  prediction  of 
maximum  responses  we  must  test  the  extension  of  the  same  general  logic  to  esti¬ 
mation  of  maximum  responses  calculated  with  the  three  different  aerodynamic 
representations.  The  results  of  this  test  are  shown  in  Tables  III  and  IV. 

TABLE  III 

RATIOS  OF  AVERAGE  RESPONSES.  OBSERVED  MAXIMUM  RESPONSES 
AND  EXPECTED  MAXIMUM  RESPONSES  FOR  PROFILE  NO .  1 


Flight  Avg.  Response  Observed  Max.  Expected  Max. 

Time  Ratios _  _ Ratios  Ratios 


Band 

II  a/ 

PP  a/ 

II 

PP 

II 

PP 

(sec.) 

Response 

PWLG 

PWLG 

PWLG 

PWLG 

PWLG 

PWLG 

1st  Bending 

50 

0.9967 

1.0039 

1.0092 

1.0034 

1.0254 

1.0036 

60 

0.9972 

1.0028 

0.9973 

1.0028 

1.0122 

1.0025 

70 

1.0024 

0.9941 

1.0129 

0.9954 

1.0110 

0.9952 

80 

0.9924 

0.9959 

1.0082 

0.9972 

1.0097 

0.9989 

90 

2nd  Bending 

0.9746 

1.0020 

0.9932 

1.0032 

1.0054 

1.0028 

50 

0.9970 

1.0049 

1.0045 

1.0043 

1.0362 

1.0046 

60 

0.9917 

1.0113 

0.9907 

1.0113 

1.0039 

1.0106 

70 

1.0091 

0.9899 

1.0178 

0 . 9911 

1.0182 

0.9910 

80 

1.0113 

0.9786 

1.0229 

0.9817 

1.0236 

0.9839 

90 

1st  Sloshing 

0.9737 

0.9993 

0.9884 

0.9999 

1.0004 

1.0000 

60 

1.0028 

1.0093 

0.9921 

1.0080 

1.0480 

1.0054 

70 

1.0117 

0.9869 

1.0472 

0.9953 

1.0376 

0.9930 

80 

1.0204 

0.9506 

1.0397 

0.9926 

1.0480 

0.9938 

90 

1.0103 

1.0045 

1.0458 

1.0029 

1.0363 

1.0026 

a/  II  denotes  instantaneous  immersion;  PP  denotes  pure  penetration;  PWLG 
denotes  penetration  with  lift  growth . 

-  50  - 


TABLE  IV 


RATIOS  OF  AVERAGE  RESPONSES.  OBSERVED  MAXIMUM  RESPONSES 
AND  EXPECTED  MAXIMUM  RESPONSES  FOR  PROFILE  NO.  2 


Flight 

Time 

Band 
(sec . ) 

Avg .  Response 
Ratios 

Observed  Max. 

Ratios 

Expected  Max. 
Ratios 

II 

PP 

II 

PP 

II 

PP 

Response 

1st  Bending 

PULG 

PWLG 

PWLG 

PWLG 

PWLG 

PWLG 

50 

1.0065 

1.0033 

1.0669 

1.0030 

1.0578 

1.0029 

60 

0.9777 

1.0033 

1.0178 

1.0024 

1.0286 

1.0025 

70 

1.0035 

0.9945 

1.0329 

0.9977 

1.0283 

0.9978 

80 

1.0015 

0.9965 

1.0140 

0.9994 

1.0264 

1.0017 

90 

2nd  Bending 

0.9928 

1.0031 

1.0048 

1.0033 

1.0114 

1.0033 

50 

1.0116 

1.0045 

1.0747 

1.0038 

1.0766 

1.0041 

60 

0.9605 

1.0134 

1.0275 

1.0104 

1.0384 

1.0107 

70 

1.0073 

0.9900 

1.0445 

0 . 9935 

1.0418 

0.9944 

80 

1.0211 

0.9811 

1.0432 

0.9846 

1.0372 

0.9889 

90 

1st  Sloshing 

0.9985 

1.0004 

1.0208 

1.0006 

1.0075 

1.0006 

60 

1.0017 

1.0075 

1.0987 

1.0002 

1.0892 

1.0008 

70 

1.0129 

0.9886 

1.0694 

0 . 9987 

1.0739 

1.0004 

80 

1.0201 

0.9577 

1.0501 

0.9972 

1.0528 

0.9990 

90 

0.9314 

1.0001 

1.0179 

1.0018 

1.0181 

1.0020 

Here,  the  penetration  with  lift  growth  results  are  used  as  a  reference  since 
this  representation  is  the  most  accurate. 


The  expected  values  of  maximum  response  used  in  the  ratios  of  Tables 
III  and  IV  were  formed  as  the  local  average  plus  one-half  the  expected  maximum 
excursion.  Where  the  expected  maximum  excursion  was  the  value  indicated  by  the 
curve  fits  in  Figs.  42  -  89.  Comparison  of  the  tabulated  ratios  for  profile 
No.  1  indicates  some  correlation  between  observed  and  expected  values  es¬ 
pecially  for  the  sloshing  response.  However,  a  significant  part  of  the  devia¬ 
tions  from  1.0  are  due  to  the  averages  which  constitute  a  sizeable  part  of  the 
response . 

The  tabulated  ratios  for  profile  No.  2,  a  weak  profile,  show  excel¬ 
lent  agreement  between  expected  and  observed  ratios. 

It  appears  that  the  differences  and  similarities  observed  in  the 
indicial  and  impulsive  responses  can  be  used  with  partial  success  to  predict 
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the  effects  of  penetration  an$  lift  growth.  The  difference  in  impulsive  excur¬ 
sions  appears  in  wind-induced  responses  as  a  change  in  extreme  excursions. 

Where  these  excursions  play  the  important  role  in  the  maximum  response  (a  weak 
wind  profile)  the  maximum  responses  are  affected  (by  penetration  primarily) 
and  to  about  one -half  the  extent  indicated  by  the  comparison  of  indicial  and 
impulsive  responses. 

In  a  moderately  severe  wind  profile  the  differences  in  averaged  re¬ 
sponses  are  likely  to  be  as  important  as  the  oscillatory  excursions.  The  dif¬ 
ferences  in  averaged  responses  are  difficult  to  predict  from  the  comparisons 
of  indicial  and  impulsive  responses.  In  addition,  it  appears  likely  that  some 
intermediate  frequencies  may  play  a  role  of  equal  significance  and  be  difficult 
to  detect  in  the  comparison  of  indicial  and  impulsive  responses. 

For  responses  which  would  be  significantly  affected  by  penetration  and 
lift  growth  it  is  likely  that  this  fact  would  be  apparent  in  the  comparison  of 
indicial  and  impulsive  responses  based  on  the  three  aerodynamic  representations. 

There  is  a  simple  and  appealing  idea  which  in  the  past  has  been  em¬ 
ployed  to  speculate  about  the  effects  of  penetration.  This  idea  is  presented 
here  and  is  shown  by  comparison  with  calculated  results  to  be  insufficient  for 
predicting  the  effects  of  penetration. 

The  displacements  in  the  bending  modes  change  sign  along  the  length 
of  the  missile.  The  generalized  forcing  function  which  drives  one  of  these 
bending  modes  takes  on  different  characteristics  when  penetration  is  neglected 
and  included. 

When  penetration  is  neglected  each  station  of  the  missile  is  immersed 
in  the  same  wind- induced  crossflow.*  The  resulting  generalized  force  for  the 
bending  mode  is  an  algebraic  sum  of  local  contributions  where  the  changes  in 
mode  shape  sign  lead  to  cancellations. 


*  The  discussion  here  pertains  only  to  wind- induced  crossflows.  The  crossflows 
due  to  the  local  transverse  body  velocities  will  always  provide  damping 
in  the  analyses  reported  here.  If  growth  of  lift  delays  were  added  to  the 
forces  from  these  crossflows  the  damping  might  be  reduced  or  eliminated. 

The  crossflows  due  to  the  local  body  angle  of  attack  are  90°  out  of  phase 
with  body  velocity  and  in  the  long  run  neither  add  nor  remove  energy  from 
the  bending  mode.  Adding  growth  of  lift  delays  to  the  forces  from  these 
latter  crossflows  could  provide  either  damping  or  undamping. 
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When  penetration  is  included  the  situation  is  best  illustrated  by  a 
unit  impulse  gust.  The  gust  crossflow  is  applied  to  successive  stations  along 
the  missile  with  delays  appropriate  for  the  time  required  to  penetrate  the 
gust.  It  is  apparent  here  that  the  energy  first  added  to  the  mode  at  the  ini¬ 
tial  penetration  may  be  augmented  or  canceled  during  the  penetration  by  sub¬ 
sequent  stations.  The  critical  factors  are  the  phase  relationships  between 
the  modal  response  and  the  succession  of  inputs.  The  character  of  previously 
calculated  responses  indicated  that  the  critical  phase  relationships  might  be 
estimated. 


Each  indicial  and  impulse  response  has  always  been  dominated  by  a 
frequency  of  the  coupled  system  which  lies  fairly  close  to  the  uncoupled  modal 
frequency.  This  led  to  the  idea  that  penetration  effects  could  be  estimated 
by  assuming  that  each  bending  mode  responds  primarily  as  an  uncoupled  mode  and 
at  its  natural  frequency.  The  comparison  of  modal  responses  with  and  without 
penetration  would  reduce  to  a  comparison  of  the  following  forms. 


R,,  (t)  pt 

_ =  /  Q„  (t)  sin  urn  (t-r)  dT  (1) 

C  J7= o  %  ^ 

Mm 


and 


6(0)  sin  ouu  (t-T)  dT 

Mm 


GL  (t)  dT 

'•m 


sin  yjT|  t 
Mm 


(2) 


where 


th 


R-  (t)  =  The  approximated 

m  eluding  penetration 


bending  mode  impulse  response  in- 


r  (t)  =  The  approximated  response  neglecting  penetration. 
Mm 

C  =  An  unknown  admittance  amplitude  coefficient. 
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QjjJt)  =  The  growth  function  for  the  m ^  bending  mode.  (Similar 
to  those  shown  in  Figs.  5  and  6.) 


=  The  uncoupled  natural  frequency. 


6(0)  =  The  unit  impulse  imposed  at  time  zero. 


sin 


=  The  assumed  modal  response  to  a  unit  impulse  imposed  at 
time  zero. 


The  approximations  R^(t)  and  have  been  evaluated  for  the 

fourth  and  second  bending  modes. m  They  are  shown  in  Figs.  36  and  37.  The  ap¬ 
proximations  for  the  second  bending  mode  response  may  be  compared  with  the  ac 
tual,  coupled  system,  impulsive  responses  shown  in  Fig.  38. 


For  both  the  fourth  and  second  bending  modes  the  approximated  re¬ 
sponses  indicate  that  incorporation  of  penetration  results  in  increased  response 
after  penetration  is  completed.  During  penetration  the  expected  sequence  of 
augmenting  and  canceling  effects  are  observed. 


The  approximated  second  mode  responses  do  not  correspond  to  the  cou¬ 
pled  system  responses.  Further,  the  major  implication  of  the  approximated  re¬ 
sponses  (increased  response  with  penetration)  is  refuted  by  the  responses  for 
the  actual  coupled  system.  The  main  reason  for  this  disagreement  appears  to 
lie  in  strong  coupling  effects,  probably  with  the  control  system  and  swivel 
engines.  These  coupling  effects  raise  the  predominant  frequency  of  response 
and  provide  a  significant  response  to  the  high  frequencies  generated  during 
penetration.  If  the  basic  idea  of  the  approximation  is  to  be  used  for  the  pre¬ 
diction  of  penetration  effects  it  will  be  necessary  to  include  some  important 
system  coupling. 


VI.  CONCLUSIONS  AND  REC CMMENDATI ON S 


Specific  conclusions  about  the  importance  of  penetration  and  lift 
growth  must  be  restricted  to  the  Saturn  C-5  without  fins  and  to  the  modes  in¬ 
cluded  in  the  analysis.  The  specific  conclusions  are: 

1.  The  calculated  responses  to  winds  are  changed  by  a  detectable  but 
insignificant  amount  when  the  aerodynamics  are  revised  to  include  penetration 
and  lift  growth. 

2.  The  change  in  responses  is  due  primarily  to  penetration;  the  addi¬ 
tion  of  lift  growth  has  very  little  effect. 
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Approximate  4th  Bending  Impulsive  Response 


■ - ii 


Fig.  36  -  Approximate  4th  Bending  Impulsive  Response  vs. 

Response  Time  -  60  sec .  F .T . 
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Fig.  37  -  Approximate  2nd  Bending  Impulsive  Response  vs 

Response  Time  - 

60  sec .  F .T. 

Actual  2nd  Bendi: 


Fig.  38  -  Actual  2nd  Bending  Impulsive  Response  vs. 

Response  Time  -  60  sec.  F.T. 
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3.  The  maximum  responses  calculated  with  instantaneous  immersion 
aerodynamics  are  conservative  (large). 

General  conclusions  are : 

1.  The  intuitive  idea  that  penetration  effects  can  be  predicted  from 
the  uncoupled  bending  mode  periods  and  corresponding  penetration  delays  is  in¬ 
correct  . 


2.  The  calculations  performed  here  provide  no  example  of  significant 
effects.  However,  it  appears  likely  that  the  cases  in  which  penetration  (or 
penetration  and  lift  growth)  plays  an  important  role  can  be  detected  by  compari¬ 
sons  of  the  indicial  and  impulsive  responses  using  aerodynamics  with  and  with¬ 
out  penetration. 

It  is  recommended  that: 

1.  The  comparison  of  indicial  and  impulsive  responses  with  and  with¬ 
out  penetration  be  used  as  a  measure  of  the  adequacy  of  instantaneous  immer¬ 
sion  aerodynamics . 

2.  Comparisons  (l  above)  be  carried  out  for  the  third  and  fourth 
bending  modes  and  for  a  Saturn  C-5  model  including  fins.  (The  fins  can  be 
simulated  by  a  conic  section  which  provides  equivalent  normal  forces  in  the 
steady  state.) 


3.  The  conclusions  of  this  report  be  checked  with  a  larger  amount  of 
wind  data  by  using  one  or  two  of  the  existing  impulsive  responses  to  calculate 
wind  responses  for  a  number  of  profiles  in  the  high  q  altitude  band. 

4.  In  the  event  a  simple  predictive  technique  for  penetration  and 
lift  growth  effects  is  sought,  consideration  should  be  given  to  refinement  of 
the  isolated  mode  idea.  The  refinement  would  include  coupling  the  mode  to  the 
important  control  frequencies. 
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APPENDIX  I 


EQUATIONS  OF  MOTION  OF  SATURN  C-5 


The  equations  of  motion  (see  ^5^)  of  a  flexible  missile  system  in 
vertical  flight  are  given  in  this  section.  Ten  generalized  coordinates  are 
considered  (see  Vol.  Ill):  lateral  translation,  yQ  ’>  rotation,  j  ,  first 
bending,  T]^  ;  second  bending,  TI2  ;  third  bending  TI3  fourth  bending,  TI4  , 
two  sloshing,  5l  ancl  ?2  J  actual  ending  deflection,  pE  >  and-  COIvtrol  deflec¬ 
tion,  Bc  (see  Fig.  39  ).  The  equations  are  valid  for  a  swivel  engine  controlled 
vehicle  where  the  swiveled  engines  account  for  four-fifths  of  the  total  thrust 
force.  Slender  body  theory  is  used  to  describe  the  generalized  aerodynamic 
forces . 


The  oscillating  propellants  are  described  by  a  mechanical  analogy 
(see  Fig.  40  )•  Only  the  motion  of  the  liquid  in  the  booster  tanks  is  investi¬ 
gated.  The  first  sloshing  mode  is  associated  with  the  furthest  aft  tank  (tank 

A) ,  while  the  second  sloshing  mode  is  associated  with  the  adjacent  tank  (tank 

B) . 

For  simplicity,  the  following  terms  of  the  equations  of  motion  (see 
[5])  are  neglected  on  the  basis  of  being  small  by  comparison: 

1.  The  rotation  of  the  missile  cross  sections  during  bending. 

2.  Icorr  (the  difference  of  mass  moment  of  inertia  of  the  frozen 
liquid  and  liquid  propellant  in  the  full  tanks  about  the  c.g.  of  the  missile). 

3.  Generalized  forces  due  to  the  flowing  propellants.  (This  elimi¬ 
nates  the  terms  containing  the  time  derivatives  of  the  mass  of  the  propellant.) 


Assuming  that  the  missile  and  atmospheric  parameters  are  constant  in 
predetermined  time  or  altitude  intervals,  the  equations  of  motion  become  for 
translation,* 


*  The  numerical  constant,  4/5,  appearing  in  (i-l),  (l-2),  and  (l-3)  as  a 

multiplier  of  0E  can  be  generalized  to  account  for  any  percentage  of 
the  total  number  of  engines  which  are  swiveled. 
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Fig.  40  -  Mechanical  Analogy  of  Sloshing  Fluids 
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i=l  l  i 

"  mSlsg^l  -  niS2Sg52  +  [^£gsE  +  ^  (xE_xg)^  ^E  =  >  (-1- 


for  mth  'bending, 
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QGm  t: 


ffiEYm(xE)  +  SEYm(xE)  +  — -  Yo 

U2  J 


QG  •* 

lE(xE-Xg)  +  SeJYjtiCxe)  +  [©E  +  SE(xE"xg)]  Ym(xE)  +  '0 


+  TA  + 


4  f 

Z  ‘  (seyiii(xe)  +  eEYm(xE^]  Yi^xE^  +  [mEYm(xE)  +  SEYm^iCE Yi^xE^ 

i=l  l 


+  — '  %  +  msiYm(xSl^l  +  mS2Ym^xS2^2 
- 

"  [SEYm(xE)  +  ®EYm(xE^J '^E  +  yo  "  [^m+Gm]^ 


+  V^Bm\  -  E  ^  % 

i=l  u 


ggjmEYm(xE)  +  ^gY^x-g)  J  +  ’  0  +  ^mwBm\i 

h  -< 

Z  |%(Wm<*E>  *  mBri(^E)Yi(XE) 


-  <3»m  ’l 


+  mslggYI|1(xsl)51  +  ms2®gYm^xS2^2 
-  [sg^x*.)  +  |  ^(xj;)]  %  =  \W 


(1-3) 


(m  =  1,2, 3, 4)  , 


for  first  sloshing, 


ft 

y0  -  (xsi‘xg)'^  +  Z  Yi(xsi)\  +  §i  +  ^sissi^i  -  Sg^ 


Z  ggyi(xSl^i  +  “ll5!  =  0  ' 


(1-4) 
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for  second  sloshing, 


4 

y0  -  (xS2"xg)0  +  H  Yi ( XS2 ) %  +  ?2  +  WS2&S2^2  ■  Sg0 
i=l 

4 

+  ^  Sg^i(xS2^%  +  U)§2^2  =  ®  •>  (l“5) 

i=l 


for  the  swivel  engine, 


-  Sgy0  +  [9E  +  SE(XE-Xg)j  0  -  H  ■  eEYi(xE^  +  SEYi^xE^|  % 

4 

+  8E8E  +  CeBe  -  CEPC  +  ggSE0  -  21  £gSEYi  ( xe ) 

i=l 

+  [kg  +  SgSE  J  3e  -  kj^C  =  0  (I_6) 

The  last  equation  of  motion  describes  the  control  and  actuator  system 
of  the  missile.  The  control  system  considered  in  this  report  utilizes  both  an 
attitude  reference  control  and  a  flow  direction  indicator  (see  QjJ ) •  The  dif¬ 
ferential  equation  describing  the  relationship  between  the  control  deflection, 
Bc  ,  the  indicated  attitude,  0i  ,  and  the  indicated  angle  of  attack,  tx±  ,  is 
given  by 


V1  V  IV  _ 

a6^C  +  a5^C  +  a4^C  +  a3®C  +  a2®C  +  al^C 

+  a0Pc  =  ao0i  +  a^i  +  b^  (l-7) 

where 

4 

0i  -  0  -  z  Yi(x0)\  (I-®) 

i=l 
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a.  =  .  yo  +  (*v-*g)  k  _  y  Yi(xv)  7)  +  0 

1  u  u  ~r,  u 


i=l 


“  Y|(xv)T)i  + 


i=l 


(1-9) 


and 


a 


v 


(1-10) 


The  quantities  afs  in  (l-7)  are  time -independent  coefficients  while 
the  gain  values  a0  }  aj_  and  b0  are  time -dependent  variables. 

The  terms  N(t)  ,  M(t)  ,  Qjj^Ct)  ,  Qj^t)  ,  Orf^t)  and  Q-^(t)  are 

generalized  aerodynamic  forces  which  are  functions  of  the  cross  flow  gust 
velocity,  vy  (see  Appendix  II). 

It  was  convenient  for  the  investigation  of  the  effects  and  importance 
of  penetration  and  growth  of  lift  on  missile  response  to  simplify  the  above 
system  of  equations  in  anticipation  of  generating  a  voluminous  amount  of  numeri¬ 
cal  results.  Thus,  for  the  study  presented  in  this  report  the  following  addi¬ 
tional  conditions  were  imposed  on  (i-l)  through  (l-9)  : 

1.  The  third  and  fourth  bending  mode  contributions  were  neglected. 

2.  The  compliance  of  the  swivel  engine  (difference  between  the 

actual  deflection  angle,  ,  and  the  control  signal,  Bq)  was  assumed  to  be 

zero.  This,  in  essence,  is  similar  to  making  the  linkage  connection  between 
the  two  variables  infinitely  rigid. 

3.  The  mass  of  the  swivel  engine  and  the  mass  moment  of  inertia  of 
the  swivel  engine  about  its  swivel  point  were  neglected. 

Incorporating  the  above  conditions  into  the  above  equations  yielded 
for  translation. 
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m  + 


+ 


+ 


+ 


u2  J  u2  u2 

mSi§i  +  ms2§2  +  Yo  “  ~  [Fl+Jo]  ^ 

|Q[Di+Hi]  T|i  +  [d2+H2]t|2  "  [f+2QFo]0 

[iri(xE)  +  QeJti-l  +  [fy2(xe)  +  0E2]ti2  -  I  FPC 


N(t)  , 


for  rotation, 


QJj. 

U2 


yo 


Icg  + 


QJ2 


-  ^  Ti!  -  ^  112 


U* 


u2 


2QF. 


-  msi(xsi-xg)§i  -  ms2(xs2-xg)§2  -  Yo  +  ^ 


[jl+F2]0 


-  [Di+Hijfl!  -  £p  [D2+H2JTI2  +  2QFi0 

-  Jpj(xE-xg)Yi(xE)  -  Yi(xe)J  -  QfiJ  % 

~  |f j(xE-xg)Y2(xE)  "  Y2^xE^J  "  ^2  j  "  mSlsg5l  '  mS2®g?2 

+  |  (xE-xg)F3c  =  M(t)  , 


for  first  bending, 


(1-11) 


(1-12) 


67 


for  second  bending, 

QGo  ,,  QGo  j  QCpi  ••  QCpp 

-r  -  -7 1*1  *  Hr  % +  Ta +  -f  ”2 

u2  u2  u2  L  u2  J 

+  mslY2(xsl)g1  +  mS2Y2(xs2)?2  +  ~~  Y0  "  jp  [52+G2J  0 

-  ^21  ^  +  ^T2t'J-g2g£2  -  •  ^2  "  2QD20  -  QA23_’ni 

+  |^T2^B2  “  QA22J Tl2  +  mSlSgY2(xSl)?l  +  mS2ggY2(xS2)§2 

-  i  FY2(xE)pc  =  Qj]2(t)  ,  (1-14) 

for  first  sloshing, 

Yo  -  (xSl‘xg)0  +  Yl(xSl)^l  +  Y2(xsi)^2  +  §1 

+  “SlSsi^l  ■  6g0  +  %Yi(xSl^l  +  %Y2(xS1^2 
+  (Ug]_5^  —  0  )  (l”l^) 
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for  second  sloshing, 


y0  ~  (xS2"xg)$  +  Yl(xS2)%  +  Y2(XS2^2  +  ?2  + 

-  Sg0  +  ggYiC^)^!  +  §gY2^xS2^2  +  tu§2l»2  =  0  * 


(1-16) 


Substituting  (l-8),  (l-9)  and  (1-10)  into  (1-7)  and  reducing  to  a  set 
of  second  order  equations  gives  for  the  control  system, 


a6P  +  ~  to  - 


al 


b0(xv-xg) 


alYi(x0)  + 


U 

t>oYl(*v) 

U 


a 


Hi  + 


alY2 


j(x0) 


fcoY2(xv) 


U 


% 


+  aipc  +  a3N  +  a5P  -  (ao+bo)0 

+  [aoYl(x0)  +  +  [aoY2(x0)  +  hoY2^]\ 


+  a0|3c  +  a2N  +  a4P  = 


(1-17) 


where 


N  -  Pc  =  0 


(1-18) 


and 


P  -  N  =  0 


(1-19) 


The  preceding  equations  of  motion  MI-11)  through  (l-19)y  describe 
the  missile  system  considered  in  the  numerical  investigation  of  the  effects  and 
importance  of  penetration  and  growth  of  lift  on  missile  response  (see  Vol.  III). 


69  - 


APPENDIX  II 


DEVELOPMENT  OF  TRANSIENT,  QUASI -STEADY  AND  STEADY  GENERALIZED 
AERODYNAMIC  FORCE  EXPRESSIONS  RESULTING  FROM  A  UNIT  STEP 
AND  UNIT  IMPULSE  WIND  PROFILE 


The  development  of  the  transient ,  quasi -steady  and  steady  generalized 
aerodynamic  force  expressions  resulting  from  a  unit  step  and  impulse  wind  pro¬ 
file  is  presented  in  this  Appendix.  The  development  is  based  on  slender  body 
theory.  For  simplicity,  the  details  of  the  analysis  are  omitted  and  the  reader 
is  referred  to  the  original  work  of  Miles  [2]  and  the  extension  to  multi -staged 
vehicles  by  Yates  [lj  for  additional  information. 

Transient  aerodynamic  force  expressions  corresponding  to  rigid  body 
and  bending  modes  of  vibration  are  presented  first  for  a  general  wind  profile. 
The  indicia!*  and  impulsive**  transient,  quasi -steady  and  steady  force  expres¬ 
sions  are  then  derived.  Quadratic  polynomials  are  used  to  curve  fit  segments 
of  the  mode  shapes  in  evaluating  the  generalized  forces  associated  with  the 
bending  modes. 

The  wind  induced  forcing  functions  which  are  compared  in  this  report 
have  both  a  geometric  and  an  aerodynamic  aspect. 

In  the  geometric  consideration  two  cases  are  used.  In  the  simplest 
case,  instantaneous  immersion,  all  stations  along  the  missile  are  assumed  to 
be  immersed  in  the  same  wind-induced  crossflow,  namely  the  wind  crossflow 
occurring  at  the  nose.  The  more  accurate  geometric  representation,  called 
penetration,  assigns  to  each  station  along  the  missile  the  wind  crossflow 
which  exists  at  the  altitude  occupied  by  the  station.  With  penetration,  the 
missile  nose  enters  a  side  gust  first  and  in  subsequent  time  successive  sta¬ 
tions  along  the  missile  length  move  into  the  crossflow. 

Two  representations  of  the  aerodynamics  are  used;  they  are  quasi¬ 
steady  and  transient.  In  the  quasi -steady  representation  the  air  forces  at  a 
missile  station  are  those  which  would  exist  if  the  local  crossflow  persisted 
unchanged  for  an  extended  time.  The  transient  representation  is  based  on 
transient  slender  body  theory  and  includes  the  growth  of  lift  with  time. 


*  Response  to  a  unit  step  wind  profile. 

**  Response  to  a  unit  impulse  wind  profile. 
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Three  types  of  wind-induced,  forcing  functions  are  assembled  using 
combinations  of  the  geometric  and  aerodynamic  representations.  The  simplest 
type  is  called  instantaneous -immersion  and  uses  the  instantaneous-immersion 
geometric  representation  with  quasi -steady  aerodynamics.  A  more  accurate  type, 
called  pure -penetration,  uses  penetration  geometries  and  quasi-steady  aero¬ 
dynamics.  The  most  accurate  forcing  functions  are  called  penetration-with-lift- 
growth.  These  functions  use  the  penetration  geometries  with  transient  aero¬ 
dynamics.* 


The  simpler  function  types,  instantaneous-immersion  and  pure-penetra- 
tion,  can  be  obtained  from  penetration -with- lift -growth  which  is  derived  first. 


I.  TRANSIENT  AERODYNAMIC  FORCES  WITH  PENETRATION  AND  LIFT  GROWTH 
FROM  A  GENERAL  CROSSFLOW  VELOCITY 


Consider  the  multi-staged  pointed  body  of  revolution  as  shown  in 
Fig.  41  .  The  cartesian  orthogonal  coordinate  system  (x-y-z)  has  its  origin  at 
the  nose.  The  vehicle  is  considered  to  be  traveling  in  the  negative  x  direc¬ 
tion  with  a  constant  velocity,  U  .  At  time  zero,  the  nose  encounters  a  side 
wind  of  magnitude  v(x,t)  directed  along  the  positive  z  axis.  Now,  from  LlJ 
and  DO  the  transient  aerodynamic  forces  corresponding  to  translational,  rota¬ 
tional  and  bending  coordinates  are  given, respectively,  by 


N(t)  =  q 


+ 


t- 

2s'(x)  r 

MR(x)  JQ 


‘±±1) 
_  MR(x) 


v(r)dT 


}  dx  , 


(II-1) 


*  The  crossflows  induced  by  missile  motions  are  in  all  cases  treated  with 

quasi-steady  aerodynamics.  These  crossflows  are  small  in  comparison  to 
the  wind-induced  crossflows  and  appear  in  the  left  hand  side  of  the 
equations  of  motion. 


71  - 


and 


+  2S'(x) 

MR(x) 


x 


u 

r u(t-  n  -?)1 

1 

f 

V  u  / 

v(T)dT  - 

0 

MR(x) 

(II-2) 


9rta(t) 


=  q. 


+  2S'(x)  P 

MR(x)  J 


‘-5 


u 


<?  5  -7) 


°L  MR(x) 


v( T )dT 


dx 


(II-3) 


where  the  notation  means  that  the  x  integration  is  performed  over 

(0,  Ut)  or  (0,L)  as  the  body  has  partially  or  totally  penetrated  the  side  wind 

profile  v^t-  •  Here  Xg  is  the  distance  from  the  nose  to  the  center  of 

gravity  and  the  sign  of  the  moment  equation  is  chosen  so  that  positive  M(t) 
gives  rise  to  a  clockwise  rotation  of  the  vehicle  as  viewed  in  Fig.  41  .  In 
(il-l),  (II-2)  and  (II-3),  primes  denote  total  differentiation  of  the  function 
with  respect  to  its  argument. 

The  exact  expression  for  f Q( t )  (see  Ql]  )  is  rather  cumbersome  to 
work  with  in  obtaining  numerical  results.  A  numerical  evaluation  of  this 
function  is  given  by  various  authors  as  are  their  expansions  for  large  and 
small  arguments  (see  [2]  and  jjEfj).  A  convenient  approximation  of  this  func¬ 
tion,  valid  for  all  values  of  the  argument,  is  given  by  Luke  W.  The  approxi¬ 
mate  mathematical  form  of  f0(T)  is 
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-0  2 1005 t  -t  -1.3063t 

f o(t)  =  1  -  0 . 0405e  -  2.7077e  +  0.0016e 


-  0.000lTe'°-21005T  -  0.3920Te"T 


-  re 


-0.61682T  (1>0204  cos  0.40731t  -  0.6574  sin  0.4073t) 


+  6-0.61682t  ^.2466  cos  0.4073r  +  0.2066  sin  0.4073t)  .  (H-4) 


This  representation  is  used  in  the  numerics  and  computer  program  presented  in 
this  report. 


XI.  TRANSIENT  AERODYNAMIC  FORCES  FOR  SATURN  C-5 
CONFIGURATION*  ENCOUNTERING  A  UNIT  STEP 
AND  UNIT  IMPULSE  WIND  PROFIIE 


A.  Unit  Step  and  Unit  Impulse  Input 

For  the  special  case  of  a  unit  step  and  unit  impulse  wind  profile, 
the  generalized  forces  corresponding  to  translational,  rotational  and  bending 
coordinates  reduce  to  a  relatively  simple  form.  We  will  first  consider  the 
indicial  transient  aerodynamic  forces. 

Let 

v(t)  =  I(t)  ,  (II-5) 

where  I  denotes  the  Heaviside  step  function.  Substituting  (11-5)  into  (il-l) 
(ll_2)  and  (II-3),  the  indicial  transient  normal  force,  moment  and  bending 
moment  become 


*  The  Saturn  C-5  configuration  considered  in  this  report  does  not  include  fins 
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Ns(t)  = 


U 


o 

U(T-  $)) 

P  S' (x)f0 ' 

\  U/ 

14R(x)  _ 

dX  ; 


Ms(t) 


(D 

-  /*  (x-xG)S'(x)f( 

u  y  b 


U 


fcDU. 


MR(x) 


(II-6) 


(II-7) 


and 


cp 

Q-p  (t)s  =  2a  f  Ym(x)S'(x)f 

■m  u  y 


U 


±1) 

°  1  MR(x) 


dx  , 


(II -8) 


where  the  subscript  s  denotes  that  these  entities  are  due  to  a  unit  step 
wind  profile. 

If  we  consider  the  unit  impulse  vind  profile 


v(t)  =  5(t)  , 


(II-9) 


where  5(t)  is  the  Dirac  delta  function,  the  impulsive  transient  normal  force, 
moment  and  bending  moment  become 


fu(T-  2  )1 

V  u  J 

-  dx  +  S'(Ut)  . 

.  MR(x) 

(II-IO) 


where 


S'(Ut)  =  S'(x) 


x=Ut 


for  0  <•  Ut  ^  L 


=  0  for  Ut  >  L 
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where 


,  (l) 


+  (UT-Xg)S'(UT) 


H(x) 


U 


hi) 


MR(x) 


dx 


(Ut-x^)S 1 (Ut )  =  (x-xg)S 1 (x) 


for  0  ^  Ut  <  L 


x=Ut 


(11-11) 


and 


=  0  for  Ut  >  L 


where 


Q71  (T)i 
*m  x 


-  q< 


(D  ,  , 

§  f  Wrs'M*° 


Kt-  5)1 


MR(x)  J 


dx 


+  S'(Ur)Ym(UT)i 


S’(UT)Ym(UT)  =  S'(x)Ym(x) 


for  0  <  Ut  <  L 


:=Ut 


(11-12) 


=  0  for  Ut  >  L 


and  the  subscript  i  denotes  that  these  entities  are  due  to  a  unit  impulse 
land  profile.  Equations  (II-6)  through  (II-8)  and  (11-10)  through  (II-12) 
define  the  indicial  and  impulsive  transient  aerodynamic  forces.  These  func¬ 
tions  will  be  used  in  the  following  section  to  describe  the  forces  for  the 
Saturn  C-5  configurations. 

Before  proceeding,  it  should  be  pointed  out  that  the  impulsive  forces 
can  be  obtained  from  the  first  time  derivative  of  the  indicial  forces.  Thus, 
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in  the  following  sections  the  indicial  forces  will  be  developed  first.  The 
latter  forces  will  then  be  differential  to  obtain  the  desired  expressions  for 
the  impulsive  forces. 


B,  Development  of  Transient  Indicial  and  Impulsive  Aerodynamic  Forces  for 
The  Saturn  C-5  Configuration 


The  basic  configuration  of  the  Saturn  C-5  missile  is  given  in 
Fig.  41.  The  nose  spike  on  the  front  of  the  missile  corresponds  to  the 
escape  tower  which  is  attached  to  the  vehicle  throughout  the  boost  flight. 

The  geometry  of  the  configuration  is  defined  in  Table  V  with  conic  character¬ 
istics  defined  by  pn  =  tan  an  for  n  =  0, 1,2,3  . 

TABLE  V 

DEFINITION  OF  SATURN  C-5  BODY  GEOMETRY 


Body  Radius 

Area  Derivative 

Region 

R(x) 

S'(x) 

Applicable 

Pox 

2rrp§x 

0  <  x  £  xG 

Ro 

0 

xQ  ^  x  ^  Xl 

P1(x-x1)+R0 

2tt0i[Bi(x-Xi)+ro] 

xi  ^  x  ^  X2 

R1 

0 

xg  ^  x  ^  X3 

Mx"x3)+r1 

2ttB2  [?2  (  x"x3 )  +Rl] 

|  X3  <  X  <  X4 

r2 

0 

X4  ^  X  <  X5 

Mx“x5)+R2 

2rrP3^P3(x-x5)+R2] 

X5  s  X  £  Xg 

r3 

0 

Xg  £  x  ^  L 

Substituting  the  appropriate  definitions  of  the  body  geometry  into 

(II-6)  yields 
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Ns(t)  =  Ng1^  (t)  +  42)(t)  +  N^(t)  +  44^(t) 


(II-13) 


where  the  components  N^t)  ,  42)(t)  ,  N<3)(t)  and  ^(t)  correspond  to 
the  normal  force  growth  on  the  first,  second,  third  and  fourth  missile  regions, 
respectively  (see  Fig.  41 ) .  These  components  are  defined  as  follows :  for 
the  first  region  (escape  tower) 


(w 


Ut^ 

:o 


N^t)  =  2q  J1'  2rrB|xf0j 
U  1 


u 


(T-  5)1 


MOoX 


dx 


0  ^  T  <  CD  , 


(11-14) 


for  the  second  region  (nose  cone) 


CUT) 

42)(t)  =  fr/  2nBipi(x-xi)+Ro]fo  | 


x- 


U 


(-5) 


Mj01(x-x1)+RoJ 
X1 


dx 


U 


^  T  <  “  , 


(11-15) 


for  the  third  region 


(xD 


43)(t)  =  Sa  f  2nB2[B2(x-x3)+R1]fo 


x3 


u 


H) 


m{B2(x-x3)+R1] 


dx 


^3 

U 


<  T  <  <*>  ; 


(11-16) 


78 


and  for  the  fourth  region  (main  booster) 


44)(t)  = 


^  J  2TTB3[p3(x-x5)+R2]f0 


u 


H) 


m[B3(x-x5)+R2^ 


l  dx 


XS 

—  £  T  <  ®  . 

u 


(11-17) 


The  limits  of  integration  in  the  above  expressions  follov  from  a 
consideration  of  the  various  penetration  times  of  the  individual  regions. 

Carrying  out  the  above  integrations  and  simplifying  for  computation, 
the  transient  indicial  normal  force  is  given  by 


Ns(t)  =  N^t)  +  42)(t)  +  i45)(t)  +  44)(t) 


(11-13) 


where  the  quantity  (T)  (n  =  1,2,3, 4)  depends  on  the  value  of  the  indepen' 

dent  variable,  t  ,  and  is  defined  by  the  following  two  expressions: 


N(n)(r)  =  2q^ 
s  U 


an(T-bn) 


1 '  2fn°Gir  ‘  fn)_ 


-  c 


n 


for  dn  ^  t  ^  d 


n 


ln  -  2fnan(T-bn)2  j^G  ;  f ^  -  G  J  ^ 


for  dn  <  t  <  »  , 
n  =  1,2, 3, 4  (11-18) 


The  coefficients  in  (II-18)  are  defined  in  Table  VI  for  the  four 
regions  of  the  C-5. 
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The  function  G(yja)  in  (II-18)  is  defined  as 


G(Y;a)  =  f  1  "■  f°^y).  dy  (II-19) 

o  (l  +  ay)3 

where  f0(y)  is  given  by  (II-4)  and  a  corresponds  to  the  product  Mf3n  (see 
Table  II).  0n  is  the  tangent  of  a  particular  region  connection  angle ,  . 

Next,  we  consider  the  moment  growth  on  the  C-5.  Substituting  the 
appropriate  definitions  of  the  body  geometry  into  (II-7)  yields 


Ms(t  )  =  M^Ct)  +  M^2)(t)  +  i43)(t)  +  l44)(r) 


(11-20) 


where 


.  .  Cv  f  u(T-  HY] 

^>(T)  =  X641)(t)  -  I1  /  *-2"fiofo 


dx 


0  ^  T  <  «  , 


(II-21) 


Ms2V) 


xgNs2)(r) 


o 

Si  f  x.2rTe1[s1(X-Xl)+Bc]f, 

Xl 


U 


(-1) 


M 


[^(x-x-^+rJ 


dx 


— i  <  7  , 

U 


(11-22) 


GO 

Ms3^(t)  =  XgNg3^ (t)  -  ia  y°  x-2rr92[02(x-x3)+R1]fo  < 

u  x3 


U 


H) 


M[p2(x-x3)+Rl^ 


dx 


*3 

U 


£  T  <  «  ,  ( 11-23) 
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and 


o 

M^4)(t)  =  xg44)(r)  -  J*  x.2np3[33(x-x5)+R2]f0 


u 


H) 


M|03(x-X5)+R2j 


l  dx 


*5 

U 


£  T  <  03  , 


(11-24) 


Now,  carrying  out  the  above  integrations  and  simplifying,  the 
transient  indicial  moment  is  given  by 


Ms(t)  =  i4^(t)  +  M^2)(t)  +  M^3)(t)  +  Mg4^  (t)  (11-20) 

where  the  quantity  Ms^t)  (n  =  1,2, 3, 4)  depends  on  the  independent  variable, 
t  ,  and  is  defined  by  the  following  two  expressions : 


M^n)(T)  =  (xg-Ubn)N(n)(T) 


U 


BfsH  (T-bn)3 

3 


l-3fnHf—  i  fn) 

\  en  J 


2Ucnen 


3f 


n 


dn  ^  T  ^  dn 


kn-2fnanU(T-bn) 


3 


H 


dn  5  t  <  03  , 

n  =  1,2, 3, 4  . 


(11-25) 


The  function  H(y;a)  in  (II-25)  is  defined  as 
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(II -26) 


H(y;a)  =  J 

O 


i-f0(y) 

(i+<?y)4 


ay 


No w,  we  will  consider  the  derivation  of  the  transient  indicial  forces 
corresponding  to  the  bending  coordinates.  From  (II-8)  it  is  seen  that  a  v 
description  of  the  mode  shapes,  Ym(x),  is  needed  in  order  to  evaluate  the 
Qjj^(T)s*s  .  The  mode  shape  data  for  the  Saturn  C-5  were  available  only  in 
discrete  numerical  form.  These  data  were  inconvenient  to  use  in  this  form 
since  they  prevented  evaluation  of  (II-8)  in  a  manner  similar  to  that  used  for 
(lI-6)  and  (II-7).  However,  by  approximating  the  mode  shapes  with  polynomials, 
the  preparation  of  (II-8)  for  computation  could  be  readily  accomplished. 

From  (II-8)  it  is  easily  seen  that  the  evaluation  of  the  integral 
requires  information  about  the  mode  shape  only  over  those  sections  of  the 
missile  which  have  a  changing  radius.  Thus  the  mode  shape  approximations  need 
only  apply  over  relatively  short  streamwise  distances.  From  an  analysis  of 
the  first  four  bending  modes  of  the  C-5  configuration,  it  was  found  that  quad¬ 
ratic  polynomials  could  be  used  efficiently  to  approximate  YmW  with  an 
acceptable  degree  of  accuracy.  In  this  instance,  the  quadratic  polynomial 
mode  shapes  yielded  values  which  were  within  1  per  cent  of  the  actual  mode 
deflections. 

Now,  let  the  mode  shapes  (m  =  1,2, 3,4)  in  the  desired  region  (see 
Fig.  41  )  be  approximated  by 


and 


Ym(x)  =  A^  +  B^x  +  C^x^  in  First  Region  (11-27) 
Ym(x)  =  A2in  +  B^x  +  C2mx2  in  Second  Region  (II-28) 
Ym(x)  =  A3m  +  B3mx  +  C3mx2  in  Third  Region  (11-29) 

Ym(x)  =  A4m  +  B^x  +  C^x2  in  Fourth  Region  (II-30) 
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where  the  A's  >  B's  3  and  C's*  are  considered  constant  for  a  discrete  flight 
time  or  altitude. 

Substituting  (II-27)  through  (II-30)  and  the  appropriate  definitions 
of  the  body  geometry  into  (II-8)  yields 

%(T)a  -  +  +  ^’(Us  (II-31) 


*  The  A's  ,  B's  and  C's  for  the  first  four  bending  modes  were  computed 
between  flight  times  of  10  and  140  sec.  at  10-sec.  intervals.  However, 
values  for  only  the  first  two  bending  modes  between  30  and  100  sec . , 
inclusive,  were  used  in  the  numerical  computation  of  the  indicial  and 
impulsive  responses. 
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«£>«.  -  (A3m  +  B3mxg)43)(T)  -  %j4S)(t) 

'Ut\ 

X4/ 


+  r  2q 
+  c3m  ~ 


U 


x2  •  2ttP2  [^2  ( x~x3 )  +Rl] f  i 


U 


(T"  u) 


M  [^2  (  x-x3 ) +Rl[j 


(  dx 


^3 

U 


<;  t  <  «  ,  (11-34) 


and 


Q^(t)3  -  (f'4n  *  E4m-<g),44)<T)  -  B4mMi4>(T) 

P 

+  c4m  J  x2-2nP3[e3(x-X5)+R2jf0 


v(t-±) 

\  U/ 


M  [p3  (  X  -X5 ) +R2J 


dx 


xs 

— — •  <  T  <  ®  . 

u 


(11-35) 


Performing  the  above  integrations  and  simplifying,  the  transient 
indicial  bending  moments  are  given  by 


«1L<T>.  “  Q^>(T)s  +  ^»)<T)b  +  Q^’(T)s  +  Q^)(T) 


(11-31) 


where 
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<V)S 


^njn  "  Cnm^ubn^2  +  xg  (2iur.  +  ^nm^nj 


»6n>(T) 


[Bnm  +  2CmUb„jM<n)(T) 


^  (T-bn)‘ 


+  —  c 
u  “M 


l-4fnl(I^;fn)]-i£(^j 

V  en  /  2  \  fn/ 


dn  5  t  s  dn 


'n 


2fnanU2(T-bn)4  [1/^  ;  fn)  l(^r^  ;  fn') 

^  \  en  /  V  en  /  J 


dn  £  t  <  «  .  ( 11-36) 


The  function  l(y;a)  in  (II-36)  is  defined  as 


i(y;a)  =  f 


¥  l-fp(y) 

(l+Oy)5 


dy 


(11-37) 


The  advantages  of  using  polynomial  curve  fits  (especially  quadratic 
polynomials)  for  the  mode  shapes ,  Ym(x)  ,  can  be  easily  seen  from  (II-36). 


The  evaluation  of  the 


can  be  made  from  the  normal  force ,  moment 


and  an  additional  integral,  l(y;a)  .  If  quadratic  polynomials  had  not  pro¬ 
vided  sufficient  accuracy  in  describing  the  Vs  ,  the  number  of  integrals 
to  evaluate  would  be  equal  to  one  less  than  the  order  of  polynomial  approxi¬ 
mation.  However,  for  the  C-5  configuration  the  transient  indicial  normal 
force,  moment  and  bending  moments  can  be  expressed  in  terms  of  integrals 
similar  to  (11-19),  (11-26),  and  (II-37). 
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As  discussed  at  the  first  of  this  Appendix,  the  normal  force,  moment 
and  bending  moments  due  to  a  unit  impulse  wind  profile  are  given  by  the  first 
time  derivative  of  the  indicial  force  and  moments.  Therefore,  differentiating 
(11-13),  (11-20)  and  (II-31)  with  respect  to  time  yields 


Nj_(T )  =  —  Ns(t)  =  (t)  +  n[2^(t)  +  n|3^(t)  +  n[4^(t)  ,  (11-38) 

dT 

M^t)  =  —  Ms(t)  =  ^^(t)  +  *42^(t)  +  (t  )  +  m|4)(t)  (11-39) 


and 


(11-40) 


where 


r 


for  dn  £  t  <  00  ; 

n  =  1,2, 3, 4  , 
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(11-41) 


M{n)(x)  =  (xg-Ubn)l4n)(T) 


r 


2fnanU(T-bn)5 

en 


H' 


U 


for  dn  s  t  £  dn 
n  =  1,2,3, 4 


-  6fnanu(T-bn): 


-  2fnanU(T-bn)3f  i_  H' 
Len 


T-dn  . 


for  dn  ^  t  <  °°  , 


n  =  1,2, 3, 4  ,  (11-42) 


and 
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^\\  ~  |^nm  “  ^nm(Ubn)2  +  xg  +  ^nm^n  ‘  Ni  (T) 

[®nm  +  2CnmU^nJ  *4  ^T) 


'run  “nrrr~n 


2anU^(T-bn)  1  -  4fnI  (~ —  *  ^n 
L  \  en  > 


2 /  „.  v  \  ^ 


2fnanU(T-bn)  T,/T-dr 


+  §ac 
u  11111 


I '  f - a  ;  fn  for  dn  5  t  £  dn  , 


n  =  1,2, 3, 4 


nanU2(T-bn)3  [I  (^2.  ;  fnj  -  I  rj±l  ;  fny 


-  2fnanU2(T-bnj  r  —  I'  (-JL  ;  fn')  “  —  I'  f  3  fn 

_en  \  en  /  en  '  en  ' 


for  dn  ^  t  <  co  } 


n  =  1,2, 3, 4  . 


(11-43) 


The  functions  G’(y;a)  ,  H'(y;a)  and  I'(y;a)  are  the  respective 
derivatives  of  G(y;ct)  ,  H(y;a)  and  i(yja)  : 


o-tvad  -  hfoh1 

(l+ay)3 


(11-44) 


H'(V;<j)  » 

(l+ay) 


(11-45) 
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i'(Yja)  = 


(11-46) 


l-f0(v) 

(l+ay)5 


When  n 


omitted. 


1 


the  terms 


and 


in  (II-41)  ,  (II-42)  and  ( 11-43),  respectively,  must  be 


The  equations  (11-13),  (II-18),  (11-20),  (11-25),  (II-31),  (11-36), 
(11-38),  (11-39)  and  ( 11-40)  are  the  fundamental  expressions  required  to 
compute  the  indicial  and  impulsive  normal  force,  moment  and  bending  moment 
growth  functions  for  the  Saturn  C-5  configuration.  These  aerodynamic  forces 
and  moments  were  used  in  calculating  the  indicial  and  impulsive  responses  of 
the  C-5  missile  system  for  the  case  ’’penetration  with  lift  growth.” 


C.  Development  of  Quasi-Steady  Indicial  and  Impulsive  Aerodynamic  Forces 
for  the  Saturn  C-5  Configuration 


The  pure -penetration  forcing  functions  are  presented  in  this  section. 
These  functions  are  based  on  penetration  geometries  and  quasi-steady  aero¬ 
dynamics. 


In  penetration  geometries  each  station  along  the  missile  experiences 
the  wind  crossflow  associated  with  the  station  altitude. 

In  the  quasi-steady  aerodynamics  the  crossflow  velocity  is  assumed  to 
be  everywhere  much  smaller  than  the  local  speed  of  sound.  Thus,  the  problem 
is  one  of  solving  Laplace* s  equation  in  each  crossflow  plane.  This  approach  is 
justified  only  in  the  case  of  low  frequency  oscillation  of  the  vehicle  at 
moderate  Mach  numbers. 

We  will  consider  first  the  penetration  normal  force,  moment  and 
bending  moment  growth  functions  due  to  a  unit  step  wind  profile.  Assuming  that 
the  local  speed  of  sound  is  essentially  infinite  with  respect  to  the  crossflow 
velocity,  then 


(II-47) 
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Substituting  (11-47)  into  (II-6),  (II-7)  and  (II-8)  gives 


(?) 

Ns(t)  =  |a  /  S'(x)I^t-  dx  ,  (11-48) 

MR— >0 

o 

Ms(t)  =  ~  ~  J*  (x-Xg)S 1  (x)l ^t-  £^dx  (11-49) 

MR — >0 

and 

(ET) 

«r  Ms  ‘  ^  J  ImMS’(x)l(T- -)**■  •  (11-50) 

u  0  \  U/ 

MR— >0 


Substituting  the  appropriate  definitions  of  the  body  geometry  into 
(II-4Q) ,  (11-49)  and  ( 11-50),  recalling  (II-27)  through  (II-30)  and  integrating 
the  results  gives :  for  the  normal  force  due  to  penetration  only. 


Ns(t)  =  Rs"L)  (t  )  +  42)(t)  +  43)(t)  +  44)(r)  (11-51) 

MR— >0  MR — ^0  MR— >0  MR— >0  MR — >0 

where 

an(T-bn)2  '  cn  for  dn  *  T  s  dn  , 

n  =  1,2, 3, 4 
hn  for  djj  s  t  <  »  , 

n  =  1,2, 3, 4  ;  (11-52) 


4n)(r)  =  2^, 
U 

MR— >0 
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for  the  moment  due  to  penetration  only, 


Ms(t)  =  M^Ct)  +  *42)(t)  +  l43)(x)  +  M^4)(t) 
MR-h»0  MR — >0  MR— >0  MR— *0  MR  >0 


where 


M^^t)  =  (xg-Ubn)4n)(r) 
MR— >0  MR— >0 


u 


2anU  * 

—  (T-bn)3 


2Uc„e 


n  n 


3f, 


Ln 


for  dn  ^  t  s  dn 
n  =  1,2, 3, 4 

for  dn  £  t  <  “  , 

n  =  1,2, 3, 4  ; 


and  for  the  bending  moment  due  to  penetration  only, 


QTl>)s  =  Qtu}(t)s  +  +  4m)(T)s  +  QV(T)s 


MR— >0 


TO 

MR— >0 


TO 

MR— *0 


TO 

MR — >  0 


MR— >0 


where 


(11-53) 


(11-54) 


(11-55) 
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MR-^O 


■^nm  "  +  xg  [®nm  + 


4n)(r) 

MR — >0 


"  [Bnm+  2CnmUbn]M^n)(T) 


MR— >0 


fanU2(T-bn) 


+  §1 


U 


'run  s 


n/^fn 

Vfn  , 


'n 


for  dn  £  t  £  dn 
n  =  1,2, 3,4 

for  dn  ^  T  <  oo  , 
n  =  1,2, 3, 4 


(II -56) 


The  penetration  normal  force ,  moment  and  bending  moment  due  to  a  unit 
impulse  wind  profile  are  obtained  from  ( 11-52),  ( 11-53)  and  (11-55),  respec¬ 
tively,  upon  differentiating  with  respect  to  time.  Thus,  for  the  normal  force 


N^t)  =  N^(t)  +  n{2^(t)  +  Np^(r)  +  n£4^(t) 
MR — >0  MR — >0  MR — >0  MR — >0  MR — >0 


where 


»'n)(-)  ■  a 

1  U 

MR— *0 


2a  (r-bn)  for 

0  for 


dn  5  r  <  dn  , 
n  =  1,2, 3, 4 

dn  £  T  <  ®  , 

n  =  1,2, 3, 4  ; 


(II-57) 


(11-58) 


for  the  moment, 


M^t)  =  41}(t)  +  m}2)(t)  +  Mp^-r)  +  m[4)(t)  (II-59) 

MR— >0  MR— ^0  MR — >0  MR — *0  MR— >0 
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where 


M{n)(T)  =  (xg-Ubn)4n)(T) 
MR — >0  MR— >0 


2anU(T-bn)^  for  dn  <.  t  £  dn  , 

n  =  1,2, 3, 4 

U  0  for  dn  s  t  <  oo  , 

n  =  1,2, 3, 4  ; 


and  for  the  bending  moments, 


%(T)i  =  +  + 

MR  — >  0  MR—>0  MR— >0 


Q\)(T)i  + 

MR— >0 


MR— >0 


where 


Cnm^n)  +  x< 


+  2C™Ub^ 
nm  n 


»|n)(T) 
MR— >0 


-  [Bnm  +  SCWJbn]  M^}(t) 

MR->0 


c 


U 


nm 


2anU2(T-bn)3  for 


0 


for 


dn  s  t  ^  dn  , 
n  =  1,2, 3, 4 

dn  s  t  <  «  , 
n  =  1,2, 3, 4  . 


(11-60) 


(11-61) 


(II -62) 
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The  equations  (ll-5l),  ( 11-53),  ( 11-55),  (11-57),  ( 11-59)  and  (II-61) 
are  the  fundamental  expressions  required  to  compute  the  penetration  indicial 
and  impulsive  force  and  moments  for  the  C-5  configuration.  The  coefficients 
given  in  these  equations  are  defined  in  Table  II.  These  forces  and  moments 
were  used  in  calculating  the  indicial  and  impulsive  responses  of  the  C-5 
missile  system  for  the  case  of  "pure  penetration." 


D.  Development  of  the  Steady  Indicial  and  Impulsive  Aerodynamic  Forces 
for  the  Saturn  C-5  Configuration 


This  section  presents  the  instantaneous-immersion  forcing  functions. 
These  functions  are  based  on  instantaneous-immersion  geometries  and  quasi¬ 
steady  aerodynamics.  In  this  geometric  consideration  every  station  along  the 
missile  experiences  the  same  wind-induced  crossflow,  namely  that  crossflow 
occurring  at  the  nose. 

We  will  consider  first  the  steady  normal  force,  moment  and  bending 
moments  due  to  a  unit  step  wind  profile.  Referring  to  (II-6),  (II-7),  and 
(II-8),  the  indicial  forces  resulting  from  instantaneous  immersion  are 


Nb(t) 


(11-63) 


and 


Ms(t) 


—  I(t) 

u 


L 

/  (x-xg)S'(x)dx 
o 


(11-64) 


(11-65) 


It  should  be  noted  that  the  aerodynamic  coefficients  in  the  above 
expressions  are  identical  to  the  coefficients  multiplying  the  translational 
velocity  terms  in  the  rigid  body  and  bending  equations  of  motion  (see  Appendix  I 
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(1-1),  (1-2)  and  (1-3)).  Thus,  (II-63),  (II-64)  and  (11-65) 
as 

can  be  rewritten* 

Ns(t)  =  I(t )  , 

u 

(11-66) 

Ms(t )  =  -  2QFl  I(t) 

(11-67) 

and 

VT>. '  ^ I(T) ' 

(11-68) 

respectively,  where  l(r)  is  the  Heaviside  step  function. 

The  instantaneous -immersion  normal  force,  moment  and  bending  moments 
due  to  a  unit  impulse  wind  profile  are  obtained  from  ( 11-66),  (11-67)  and 
(11-68)  upon  differentiating  with  respect  to  time: 

%(t)  =  2^°  6(t) 

(11-69) 

%(t)  =  -  2QFl  5(t) 

(11-70) 

and 

■  tt  8(t> 

(11-71) 

*  The  evaluation  of  Ns(t)  ,  Ms(t)  and  Op  (t)s  can  also  be  obtained  from 
(11-52),  (11-54)  and  (11-56)  by  setting  t  equal  to  infinity. 


where  6(t)  is  the  Dirac  delta  function. 

For  t  >0  ,  %(r)  ,  M^(t)  and  T) i  are  identically  equal  to 

zero  (see  (11-58),  (II-60)  and  (11-62)  for  t  =  »  )  . 

The  equations  (11-66)  through  (ll-7l)  are  the  fundamental  expressions 
required  to  compute  the  instantaneous-immersion  indicial  and  impulsive  force 
and  moments.  These  forces  were  used  in  calculating  the  indicial  and  impulsive 
responses  of  the  Saturn  C-5  missile  system  for  the  case  of  "instantaneous 
immersion. " 
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APPENDIX  331 


CALCULATION  OF  INITIAL  AND  STEADY -STATE  CONDITIONS  RESULTING 
FROM  A  UNIT  IMPULSE  AND  UNIT  STEP  WIND  PROFILE 


The  numerical  solution  for  the  indicial  and  impulsive  responses  of 
the  vehicle  requires  a  knowledge  of  the  initial  conditions  imposed  by  the  var¬ 
ious  forcing  functions  (see  Appendix  XI).  Since  the  Runge  Kutta  method  of 
solution  is  used,  the  initial  conditions  for  the  generalized  coordinates  and 
their  first  derivatives  are  needed. 

The  steady-state  values  for  the  indicial  and  impulsive  responses  of 
the  vehicle  are  required  for  computer  logic  which  terminates  the  integration 
when  response  is  sufficiently  close  to  its  steady  value  (see  Vol.  III). 

A  procedure  is  given  in  this  Appendix  for  computing  the  initial  and 
steady-state  conditions  for  the  impulsive  and  indicial  responses  of  the  vehicle. 
For  simplicity,  the  method  is  described  for  a  system  defined  by  the  following 
generalized  coordinates  (see  Appendix  I):  translation,  rotation,  first  and 
second  bending,  first  and  second  sloshing  and  control  deflection.  The  procedure 
is  general  in  nature,  however,  and  can  be  applied  to  a  system  with  more  degrees 
of  freedom. 

The  initial  and  steady-state  conditions  are  obtained  through  usage 
of  Laplace  transform  techniques.  The  initial  and  steady-state  conditions  are 
fouhd  in  the  limit  of  the  transform  as  the  Laplace  variable  approaches  infinity 
and  zero,  respectively.  The  above  conditions  are  derived  for  six  cases  of 
forcing  functions:  unit  impulse  and  step  considering  penetration  with  lift 
growth,  pure  penetration  and  instantaneous- immersion  effects . 


I.  CALCULATION  OF  INITIAL  CONDITIONS 

Writing  the  equations  of  motion  [see  (i-ll)  through  (1-19 )j  as  a  set 
of  first  order  equations,  we  find 

{*}-{«}  (III-1) 
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where  the  elements  and  (i,j  =  1,2,.... ,18)  of  the  square  matrices 

[X]  and  [jj] ,  respectively,  correspond  to  the  coefficients  on  the  left  hand  side 
of  (i-ll)  through  (1-19), 


(III-2) 


and  {c}  is  a  column  matrix  of  forcing  functions .  Six  different  sets  of  forcing 
functions  were  considered  in  the  numerical  investigation  presented  in  this  re¬ 
port  . 


For  conveniences  of  comparison  the  elements,  c^  ,  (i  =  1,2,...., 18) 
of  <[c}  are  given  in  Table  VH  for  a  unit  impulse  wind  profile  {vy( t )  =  6(t)j 
and  a  unit  step  wind  profile  {v^(t)  =  l(t )j»  considering  penetration  with  lift 
growth  (FWLG),  pure  penetration  (PP)  and  instantaneous -immersion  (II)  effects. 

Expressions  for  the  aerodynamic  quantities  in  Table  VII  are  given  in 
Appendix  II. 
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TABLE  VII 


ELEMENTS  OF  FORCING  FUNCTION  MATRIX 


Unit 

Impulse  Wind  Profile 

Unit  Step  Wind  Profile 

i 

FWLG 

■PP 

II 

FWLG 

PP 

II 

i 

0 

0 

0 

0 

0 

0 

2 

0 

0 

0 

0 

0 

0 

3 

0 

0 

0 

0 

0 

0 

4 

0 

0 

0 

0 

0 

0 

5 

0 

0 

0 

0 

0 

0 

6 

0 

0 

0 

0 

0 

0 

7 

0 

0 

0 

0 

0 

0 

8 

0 

0 

0 

0 

0 

0 

9 

0 

0 

0 

0 

0 

0 

10 

VT> 

N.(t) 

m->o 

Vt) 

Hs(t) 

HS(T) 

MR— >  0 

Ns(t) 

ii 

M^t) 

Ml(t) 

MR-?0 

M±(t) 

Ms(t) 

M  (t) 
sv  ' 

MR— ?  0 

Ms(t) 

12 

y  t)i 

\(T)i 
MR—?'  0 

Si('T^i 

\^s 

Yt>s 

MR— ?  0 

\(T)s 

13 

\<T)i 

MR—?  0 

VT)i 

VT)s 

VT>s 

MR— ?  0 

VT)s 

14 

0 

0 

0 

0 

0 

0 

15 

0 

0 

0 

0 

0 

0 

16 

0 

0 

0 

0 

0 

0 

17 

0 

0 

0 

0 

0 

0 

18 

Y>6(t) 

t>o6(T) 

V(t) 

Y)I(t) 

v(T) 

V(T) 

u 

u 

u 

u 

u 

u 
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Now,  the  initial  conditions,  <[q(0)}.,  can  he  found  readily  by  using 
Laplace  transform  techniques.  Taking  the  Laplace  transform  of  (III-1)  yields 


where 


and 


[ab]  {l}  -  {c} 

(III-3) 

iQ  1 

II 

H 

3 

(III-4) 

[ab]  =  s[a]  +  [b] 

(III-5) 

The  elements  c±  (i  =  1,2, ...,18)  of  the  transformed  forcing  function 
matrix  {«}  are  given  in  Table  VIII  for  the  six  different  sets  of  forcing  func¬ 
tions  .  The  symbol  &  denotes  "the  order  of  i 

The  initial  conditions ,  *  <[q(o)]>,  are  now  obtained  by  solving  (ill- 3) 
for  each  dependent  variable  by  Cramer’s  rule,  multiplying  the  ratio  of  two 
determinants  by  s  and  taking  the  limit  as  s  — >  co  . 

The  initial  conditions  imposed  by  a  unit  impulse  wind  are  given  in 
Table  IX  for  the  cases  where  penetration  with  lift  growth  and  pure  penetration 
effects  are  considered. 


*  If  f(s)  is  given  as  the  Laplace  transform  of  f(t)  ,  then  the  initial 
condition,  f(+o)  ,  is  found  from 


lim  sf(s)  =  f(+o) 

S - >  co 


where  s  is  the  Laplace  transform  variable. 
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TABLE  mi 


ELEMENTS  OF  THE  TRANSFORMED  FORCING  FUNCTION  MATRIX 


Unit  Impulse  Unit  Step 

Wind  Profile _  Wind  Profile 


-L 

FWLG  and  PP 

n 

FWIG  and  PP 

II 

i 

0 

0 

0 

0 

2 

0 

0 

0 

0 

3 

0 

0 

0 

0 

4 

0 

0 

0 

0 

5 

0 

0 

0 

0 

6 

0 

0 

0 

0 

7 

0 

0 

0 

0 

8 

0 

0 

0 

0 

9 

0 

0 

0 

0 

10 

2QFc 

u 

<30 

fQfo  1 

U  s 

11 

<30 

-2QFX 

u 

<70 

-2QFx  1 

U  s 

12 

oil ) 

2QDX 

U 

<?0 

2QDX  1 

U  s 

13 

<30 

2QD2 

U 

<0 

2QP2 

U  s 

14 

0 

0 

0 

0 

15 

0 

0 

0 

0 

16 

0 

0 

0 

0 

17 

0 

0 

0 

0 

18 

b 

b 

b  n 

o 

o 

o  1 

Do  1 

U 

U 

U  s 

U  s 
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TABLE  IX 


INITIAL  CONDITIONS  IMPOSED  BY  A  UNIT  IMPULSE  WIND 
FOR  THE  CASES  OF  PENETRATION  WITH  LIFT  GROWTH 
AND  PURE  PENETRATION 


y0(°)  -  0 

y0(o) 

4(o)  =  o 

k  o) 

£ 

o 

It 

o 

%(0) 

o 

II 

o 

V°) 

§x(o)  =  o 

§x(0) 

§2(0)  =  0 

Sg(o) 

pc(0)  =  0 

M°) 

S3 

O 

II 

O 

N(0) 

P(0)  =  0 

P(0) 

The  initial  conditions  imposed  by  a  unit  impulse  wind  for  the  case 
of  instantaneous  immersion  are: 


y0(°)  =  0 

jZS(O)  =  0 
^(0)  =  0 
T)2(o)  =  0 
gpCo)  =  o 
S2(0)  =  o 
Pc(0)  =  0 


N(0)  =  0 

P(0)  =  0 

Efc(o)  =  0 
N(0)  =  0 


(III-6) 
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•  •  •  •  • 

In  addition,  solutions  for  yo(0)  ,  0(0)  ,  Tjj^O)  ,  TjgC 0)  ,  ^(0) 
and  i2( °)  are  obtained  from  the  matrix  equation 

{z(o)}  »  [d]'1^}  (in-7) 

where 


The  elements,  ,  of  the  square  matrix  £dJ  correspond  to  the 

inertia  coefficients  in  (i-ll)  through  (1-16)  (see  Appendix  I). 

Finally,  the  initial  conditions  imposed  by  a  unit  step  wind  are  given 
in  Table  X  for  the  cases  where  penetration  with  lift  growth,  pure  penetration 
and  instantaneous- immersion  effects  are  considered. 
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TABLE  X 


INITIAL  CONDITIONS  IMPOSED  BY  A  UNIT  STEP  WIND  FOR 
THE  CASES  OF  PENETRATION  WITH  LIFT  GROWTH,  PURE 
PENETRATION  AND  INSTANTANEOUS  IMMERSION 


yQ(o)  =  o 
4(0)  =  0 
\(0)  =  0 
Vq(o)  =  o 

§3.(0)  =  o 
§2(0)  =  0 
pc(0)  =  0 
N(0)  =  0 

P(0)  =  0 


yQ(o)  =  o 
$( o)  -  0 
%(0)  =  0 
T12(0)  =  0 

§3.(0)  =  o 

§2(0)  =  0 

Sfc(o)  =  0 
N(0)  =  0 

P(0)  =  0 


II,  CALCULATION  OF  STEADY-STATE  VALUES 


The  procedure  for  computing  the  steady-state  values  of  the  impulsive 
and  indicial  responses  follows  the  analysis  given  in  Section  I  except  for  the 
following  changes: 

1.  The  first  order  system  of  equations  of  motion  [see  (ill-l)]  are 
written  in  terms  of  y0  and  yQ  instead  of  yQ  and  yQ  to  avoid  an  inde¬ 
terminant  form. 

2.  The  steady  state  values  of  the  impulsive  responses  are  the  same 
for  the  three  aerodynamic  environments  considered.*  The  same  is  true  of  the 
steady  state  values  of  the  indicial  responses.  Thus,  the  steady  state  values 
for  the  impulsive  and  indicial  responses  will  be  found  from  the  case  of  instan 
taneous  immersion. 
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3.  The  steady  state  conditions  are  obtained  when  the  Laplace  trans¬ 
form  variable  goes  to  zero  in  the  limit .** 

Omitting  the  details ,  the  steady  state  values  of  the  impulsive  and 
indicial  responses  are  given  in  Tables  VII  and  VIII,  respectively. 


TABLE  XI 

STEADY  STATE  VALUES  FOR 
IMPULSIVE  RESPONSES 


yo(a>)  =  0 

o-: 

IT 

1! 

o 

$((<»)  =  0 

0(co)  =  0 

%(“)  =  0 

%(“>)  »  0 

Tfe(»)  =  0 

T^C00)  =  0 

§!(“)  =  0 

^(co)  =  0 

§2(o>)  =  0 

52(a>)  =  0 

Pc(eo)  =  0 

j3c(a>)  =  0 

N(»)  =  0 

N(oo)  =  0 

P(»)  =  0 

P(~)  =  0 

TABLE  XII 

STEADY  STATE  VALUES  FOR 
INDICIAL  RESPONSES 


yQ(“) 

=  1 

y0(m)  =  o 

0(oo) 

=  0 

IT 

11 

o 

%(00) 

=  0 

%(“)  =  0 

=  0 

o 

11 

T 

S^00) 

=  0 

§!(“)  =  0 

?2  (“) 

=  0 

§2(oc)  =  0 

ec(°0 

=  0 

Pc(°°)  =  0 

N(eo) 

=  0 

N(co)  =  0 

p(«) 

=  0 

P(co)  =  0 

*  The  steady  state  values  of  impulsive  forcing  functions,  given  in  Table  III, 
are  the  same  for  the  three  aerodynamic  environments .  The  same  is  true  for 
the  steady  state  values  of  the  indicial  forcing  functions .  Since  the 
forcing  function  vector  is  the  only  quantity  which  changes  in  the  set  of 
equations  [see  (Ill-lj]  for  the  different  aerodynamic  environments,  it 
is  easily  seen  that  the  above  statement  is  correct. 

**  If  f(s)  is  the  Laplace  transform  of  f(t)  ,  then  the  steady  state  value, 
f(+®)  ,  is  found  from 


lim  sf ( s )  =  f(+«) 
s - >0 
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APPENDIX  IV 


PRESENTATION  OF  EXTREME  EXCURSION  AND 
AVERAGE  RESPONSE  PLOTS 


This  Appendix  contains  the  extreme  excursion  and  local  average 
response  plots  discussed  in  Section  V-  The  method  of  analysis  used  to  obtain 
these  plots  is  discussed  in  Section  IV. 
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Fig.  42  -  Extreme  Excursion  of  1st  Bending  vs.  Reduced  Cumulative 

Probability  -  50  sec.  F.T. 


Fig.  43  -  Extreme  Excursion  of  1st  Bending  vs.  Reduced  Cumulative 

Probability  -  60  sec.  F.T. 


Fig.  44  -  Extreme  Excursion  of  1st  Bending  vs.  Reduced  Cumulative 

Probability  -  70  sec.  F.T. 


Fig.  45  -  Extreme  Excursion  of  1st  Bending  vs.  Reduced  Cumulative 

Probability  -  80  sec.  F.T. 
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Fig.  46  -  Extreme  Excursion  of  1st  Bending  vs.  Reduced  Cumulative 

Probability  -  90  sec.  F.T. 


Ill 


Fig.  4?  -  Extreme  Excursion  of  2nd  Bending  vs.  Reduced  Cumulative 

Probability  -  50  sec.  F.T. 


Fig.  48  -  Extreme  Excursion  of  2nd  Bending  vs.  Reduced  Cumulative 

Probability  -  60  sec.  F.T. 
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Fig.  49  -  Extreme  Excursion  of  2nd  Bending  vs.  Reduced  Cumulative 

Probability  -  70  sec.  F.T. 


Fig.  50  -  Extreme  Excursion  of  2nd  Bending  vs.  Reduced  Cumulative 

Probability  -  80  sec.  F.T. 
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Fig.  51  -  Extreme  Excursion  of  2nd  Bending  vs.  Reduced  Cumulative 

Probability  -  90  sec.  F.T. 
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Fig.  52  -  Extreme  Excursion  of  1st  Sloshing  vs.  Reduced  Cumulative 

Probability  -  60  sec.  F.T. 


Fig.  53  -  Extreme  Excursion  of  1st  Sloshing  vs.  Reduced  Cumulative 

Probability  -  70  sec.  F.T. 
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(PWLG) 


Profile  No.  1 
80  sec.  F.T.  Band 


Fig.  54  -  Extreme  Excursion  of  1st  Sloshing  vs.  Reduced  Cumulative 

Probability  -  80  sec.  F.T. 


Fig.  55  -  Extreme  Excursion  of  1st  Sloshing  vs.  Reduced  Cumulative 

Probability  -  90  sec.  F.T. 
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Fig.  56  -  Extreme  Excursion  of  2nd  Sloshing  vs.  Reduced  Cumulative 

Probability  -  60  sec.  F.T. 


Fig.  57  -  Extreme  Excursion  of  2nd  Sloshing  vs.  Reduced  Cumulative 

Probability  -  70  sec.  F.T. 
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Fig.  59  -  Extreme  Excursion  of  2nd  Sloshing  vs.  Reduced  Cumulative 

Probability  -  90  sec.  F.T. 
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2.0 


Reduced  Cumulative  Probability 

Fig.  64  -  Extreme  Excursion  of  1st  Control  vs.  Reduced  Cumulative 

Probability  -  90  sec.  F.T. 
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Fig.  69  -  Extreme  Excursion  of  1st  Bending  vs.  Reduced  Cumulative 

Probability  -  70  sec.  F.T. 


Fig.  70  -  Extreme  Excursion  of  1st  Bending  vs.  Reduced  Cumulative 

Probability  -  80  sec.  F.T. 
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Fig.  71  -  Extreme  Excursion  of  1st  Bending  vs.  Reduced  Cumulative 

Probability  -  90  sec.  F.T* 
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Pig.  74  -  Extreme  Excursion  of  2nd  Bending  vs.  Reduced  Cumulative 

Probability  -  70  sec.  F.T. 


Fig.  75  -  Extreme  Excursion  of  2nd  Bending  vs.  Reduced  Cumulative 

Probability  -  80  sec.  F.T. 
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Reduced  Cumulative  Probability 

Fig.  76  -  Extreme  Excursion  of  2nd  Bending  vs.  Reduced  Cumulative 

Probability  -  90  sec.  F.T. 
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Reduced  Cumulative  Probability 

Fig.  77  -  Extreme  Excursion  of  1st  Sloshing  vs.  Reduced  Cumulative 

Probability  -  GO  sec.  F.T. 


Reduced  Cumulative  Probability 

Fig.  78  -  Extreme  Excursion  of  1st  Sloshing  vs.  Reduced  Cumulative 

Probability  -  70  sec.  F.T. 
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Reduced  Cumulative  Probability 


Fig.  79  -  Extreme  Excursion  of  1st  Sloshing  vs.  Reduced  Cumulative 

Probability  -  60  sec.  F.T. 


Reduced  Cumulative  Propability 


Fig.  80  -  Extreme  Excursion  of  1st  Sloshing  vs.  Reduced  Cumulative 

Probability  -  90  sec.  F.T. 
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Reduced  Cumulative  Probability 


Fig.  81  -  Extreme  Excursion  of  2nd  Sloshing  vs.  Reduced  Cumulative 

Probability  -  60  sec.  F.T. 


Reduced  Cumulative  Probability 

Fig.  82  -  Extreme  Excursion  of  2nd  Sloshing  vs.  Reduced  Cumulative 

Probability  -  70  sec.  F.T. 
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Reduced  Cumulative  Probability 

Fig.  84  -  Extreme  Excursion  of  2nd  Sloshing  vs.  Reduced  Cumulative 

Probability  -  90  sec.  F.T. 
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Reduced  Cumulative  Probability 

Fig.  88  -  Extreme  Excursion  of  1st  Control  vs.  Reduced  Cumulative 

Probability  -  80  sec.  F.T. 
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Reduced  Cumulative  Probability 

Fig.  e9  -  Extreme  Excursion  of  1st  Control  vs.  Reduced  Cumulative 

Probability  -  90  sec.  F.T. 
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Average  Response  -  1st  Bending 


Average  Respon; 


Profile  No.  1 


Average  Response  of  1st  Sloshing  vs.  Flight  Time 


Flight  Time  (sec.) 


Fig.  93  -  Average  Response  of  2nd  Sloshing  vs.  Flight  Time 


Average  Response  -  2nd  Bending  x  10  ^  Average  Response  -  1st  Bending 


Flight  Time  (sec.) 


Fig.  95  -  Average  Response  of  1st  Bending  vs.  Flight  Time 


Flight  Time  (sec.) 


Fig.  96  -  Average  Response  of  2nd  Bending  vs.  Flight  Time 


Average  Response  -  2nd  Sloshing 


Average  Response  -  1st  Control 


Flight  Time  (sec.) 


Fig.  99  -  Average  Response  of  1st  Control  vs.  Flight  Time 


NASA -Langley,  1965  CR-3  ^6 
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"The  aeronautical  and  space  activities  of  the  United  States  shall  be 
conducted  so  as  to  contribute  ...  to  the  expansion  of  human  knowl¬ 
edge  of  phenomena  in  the  atmosphere  and  space.  The  Administration 
shall  provide  for  the  widest  practicable  and  appropriate  dissemination 
of  information  concerning  its  activities  and  the  results  thereof. ,f 

— National  Aeronautics  and  Space  Act  of  1958 
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